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1 .  INTRODUCTION 


The  recent  trend  to  high  speed  aircraft  and  missiles 
has  subjected  these  vehicles  to  the  additional  phenomenon  of 
aerodynamic  heating.  This  heating  must  be  considered  in  any 
analysis  of  the  vehicle’s  structural  integrity.  Since  exper¬ 
imental  testing  of  small  scale  models  can  play  an  important 
role  in  supplementing  and  substantiating  theoretical  analysis, 
the  need  for  a  good  understanding  of  the  general  similarity 
laws  of  this  combined  aerodynamic,  thermal,  and  elasticity 
problem  is  apparent. 

The  present  report  reviews  some  of  the  basic  r.imi- 
larity  parameters  of  this  combined  aerothermoelastic  problem 
and  assesses  their  general  significance.  Also,  it  will  at¬ 
tempt  to  indicate  means  of  dealing  with  some  of  the  conflict¬ 
ing  requirements  of  the  similarity  parameters.  For  simplicity, 
the  present  study  is  limited  to  moderate  temperatures  (below 
about  1000°F  -  i«e.,  M  <  3  =  5 )  and  considers  only  elastic  be¬ 
havior  of  the  material. 

A  general  summary  of  the  overall  high  temperature 
aircraft  structural  problem  can  be  found  in  Ref.  1,  together 
with  some  comments  regarding  the  use  of  models  in  high  temper¬ 
ature  structural  design.  The  specific  prohl&n  of  aerothermo¬ 
elastic  similarity  laws  has  also  been  studied  by  several  in¬ 
vestigators  during  the  past  six  years.  Refs.  2  through  10 
give  a  list  of  some  of  the  approaches  to  this  problem. 

1.1  The  Aerothermoelastic  Model 

For  flight  at  Mach  numbers  less  than  about  1.5 
where  one  c  u  neglect  aerodynamic  heating,  the  testing  of 
cold  aeroelastic  models  (in  particular  flutter  models),  has 
emerged  as  an  indispensible  aid  in  designing  structures  to 
avoid  aeroelastic  phenomena.  There  are  two  basically  differ- 
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ent  approaches  to  the  use  of  models  in  aeroelastie  development 
work.  In  one,  the  model  tests  are  designed  to  evaluate  coef^- 
ficients  to  be  inserted  into  the  theories.  In  the  other,  the 
model  tests  are  designed  as  analogues  of  the  full  scale  prob¬ 
lem.  The  latter  approach  is  particularly  desirable  in  proof 
testing  specific  aircraft  configurations,  and  it  is  here  that 
a  thorough  understanding  of  the  similarity  laws  is  necessary. 

At  the  higher  flight  Mach  numbers,  it  becomes  neces¬ 
sary  to  include  the  effects  of  aerodynamic  heating  in  any 
aeroelastie  model  test  conducted  to  represent  the  full  scale 
aircraft.  One  means  of  doing  this  is  to  first  compute  the 
loss  in  stiffness  due  to  the  aerodynamic  heating,  then  to 
construct  the  model  to  conform  to  this  reduced  stiffness, 
and  finally  to  test  this  reduced  stiffness  model  in  a  conven¬ 
tional  low  stagnation  temperature  wind  tunnel  at  the  proper 
Mach  number.  This  approach  requires  analytical  computation 
of  heating  effects  and  also  generally  assumes  separation  into 
separate  aerothermal  and  aeroelastie  problems  with  no  feed¬ 
back  interaction  between  the  two.  See  the  discussion  in  Ref.  1. 

A  more  direct  approach  to  this  problem  would  be  to 
construct  a  model  and  insert  it  directly  into  a  high  stagna¬ 
tion  temperature  wind  tunnel.  This  aerothermoelastic  model 
is  thus  subjected  to  the  heating  and  airloads  simultaneously. 
This  approach  requires  a  thorough  knowledge  of  the  similarity 
conditions  for  aerodynamic  heating  as  well  as  aeroelasticity. 

As  will  be  seen,  the  satisfaction  of  all  the  similarity  condi¬ 
tions  for  an  aerothermoelastic  model  is  not  generally  possible 
except  for  a  scale  ratio  of  1  :  1.  Hence  either  some  simi¬ 
larity  condition  will  have  to  be  relaxed  in  extrapolating  the 
model  to  full  scale  or  the  model  test  will  have  to  be  used  in 
the  sense  of  evaluating  coefficients  t;o  be  inserted  into  vari¬ 
ous  theories. 
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In  addition  to  these  complete  aerothermoelastic 
models,  "restricted  purpose"  models  might  be  used  to  advantage 
to  investigate  certain  well-defined  phenomena.  Aerodynamic 
heat  transfer  models,  cold  aeroelastic  models,  structural 
models,  and  cold  aerodynamic  stability  derivative  models  are 
some  well  known  examples  of  such  ''restricted  purpose"  models 
that  can  yield  significant  data  on  specific  problem  areas. 
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2.  GENERAL  AEROTHERMQELAST I C  SIMILARITY  LAWS 

Various  studies  of  the  similarity  laws  associated 
with  models  placed  in  high  stagnation  wind  tunnels  are  given 
in  Refs.  2  to  10.  Generally,  these  studies  are  made  either 
by  applying  the  Buckingham  7T  theorem  of  dimensional  analysis 
as  in  Refs.  6  and  9,  or  by  examining  the  appropriate  governing 
equations  in  non-dimensional  form,  as  in  Refs.  4,  5,  and  7. 

The  latter  approach  seems  more  flexible,  particularly  since 
the  significance  of  the  various  non-dimensional  parameters 
can  be  better  ascertained.  This  approach  will  be  utilized 
in  this  section. 

Following  Ref.  5,  the  general  aerothermoelastic 
problem  can  be  divided  into  three  parts,  namely,  one  dealing 
with  the  aerodynamic  flow,  one  dealing  with  the  heat  conduc¬ 
tion  to  the  interior,  and  a  third  part  dealing  with  the  ther¬ 
mal  stresses  and  external  loads  applied  to  the  structure. 

The  similarity  relations  associated  with  each  of  these  parts 
will  be  treated  in  turn. 


2.1  Aerodynamic  Similarity 

* 

The  two-dimensional  flow  of  a  compressible  viscous 
heat  conducting  perfect  gas  is  characterized  by  the  following 
equations. 


tL 

Dt 


t  <L  (p*)  +  j L  (p+J) 

e?X  Pi 


O 


(2.1; 


■* 

p  d u  /.  u.  is  f  w  dti  -  _  df 
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Ll 

3 


t  ?u  iS  iS 

x 1 1  2%  ' 


r* 

ir  / 


+  fJL  ijl 

Pi  Lr  d* 


iS. 

pi 


(2.2) 


Two-dimensional  flow  is  assumed  here  for  simplicity  in 
writing.  The  three-dimensional  situation  will  not  intro 
duce  any  additional  parameters. 
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where  u,  w  represent  the  velocity  components,  T  the  tempera¬ 
ture,  P  the  density,  p  the  pressure,  n  the  viscosity,  k  the 
heat  conductivity,  cp  the  specific  heat  at  constant  pressure, 
R  the  universal  gas  constant,  and  m  the  molecular  weight  of 
the  gas. 

If  now  the  following  non-dimensional  parameters  are 
introduced, 

X  --  ?  ---  Z 
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where  the  subscript  o  denotes  some  reference  value  and  the 
subscript  <»  denotes  the  free  stream  value,  the  above  equations 
will  reduce  to. 


PcoV\(L__\il  +  <L(fu)  -h  a  (fwj\  *  O 

L  ^  ** 

p  [7jL  )  IE:  +u.3“  +w  3*}  =  -/ft  )  i£  +/H*  )f  Ai.fi  p 

r  vvi.hi  a*  3%  J  (fmvy  3*  ' 


(2.7) 


Ti  V  **  7% 


dJL-U^i 

&  si 


(2.8) 
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+(*.  )[£.(* il)  (k  iZ 
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(2.10) 


^£7  »'  (ft  fb  7*  ]  /  ^0  )  />  7" 
r  X^TfiTJ  U  / 


f  T 


(2.11) 


Examination  of  the  above  equations  reveals  that  if 
certain  non-dimensional  parameters  are  the  same  for  a  model 
and  a  prototype,  then  Eqs,  (2.7)  to  (2.11)  will  be  identical 
for  both  these  cases  and  will  have  identical  non-dimensional 
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solutions  for  the  unknowns  u,  w,  T,  p  and  pe  Hence  measure 
ments  of  ii,  w ,  T,  P,  and  p  on  a  model  can  then  be  used  to 
find  u,  w,  T,  P ,  and  p  on  the  prototype  by  application  of 

Eqs,  (2.6),  These  non-dimensional  parameters  that  must  be 
identical  are, 

A 

U-to 

fa>  _ 

fcu1 


/U*  _  / 

JoL 


Too 

To 

,  I 


Thus  similarity  in  the  aerodynamic  flow  and  heating  will  pre 
vail  if  the  following  nine  parameters  are  maintained  equal. 

Free  stream  Mach  number.  M  —  a & 

Free  stream  Reynolds  number,  A 

H" 


12 
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Free  stream  Prandtl  number. 

Ratio  of  Specific  heats. 

Unsteady  Flow  parameter. 

Reference  Temperature  ratio, 

_  _  af 

/d  f  }  j  •£ 

Usually  the  aerodynamic  heating  model  test  will  be 
made  in  the  same  medium  (air)  as  the  prototype*  In  this 
case,  the  y  ' s  are  the  same.  Also,  the  Fr  generally  doesn't 
vary  much  over  a  large  temperature  range.  The  unsteady  term 
^/l  and  reference  temperature  ratio  ~r'/ Tm>  sarve  to  define 

'  it 

the  reference  time  tQ  and  the  reference  temperature  T0  re¬ 
spectively.  The  quantities  Cp,  H,  k  depend  only  on  the  tem¬ 
perature  T.  If  the  free  stream  temperature  Is  chosen  the 
same  for  both  model  and  prototype  in  the  same  gas,  these  non- 
dimensional  quantities  will  always  be  the  same.  If  the  T»  's 
are  somewhat  different,  these  conditions  can  still  be  main¬ 
tained  over  a  low  range  of  temperatures  since  cp  does  not 

vary  much  between  0  -  1000°F  and  k  and  \x  can  be  expressed  as 

o(  9 

k  ^T  and  g.-T  over  this  limited  range. 

Under  all  these  above  conditions  the  aerodynamic 
heating  similarity  reduces  to  the  requirements  mainly  of 
equal  Mach  number  M  and  equal  Reynolds  number  Re  for  the 
flow. 

In  addition  to  the  differential  equations,  the 
boundary  conditions  should  also  be  investigated.  The  bound¬ 
ary  conditions  at  the  surface  of  the  body  are. 


It  will  be  convenient  when  later  dealing  with  heat  conduction 
similarity  to  choose  some  reference  temperature  other  than 
the  free  stream  temperature. 
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(2 .  14 ) 


where  q ^  is  the  heat  flux  per  unit  area  and  n  is  the  normal 
direction  to  the  body  surface.  Introducing  the  non-dimen¬ 
sional  quantities  of  Eq.  (2,6)  into  the  above  will  result  in 
the  requirement  that  prototype  and  model  have  the  same  exter- 

•fa 

nal  geometric  boundary  shape,  and  that  the  air  temperature 
be  equal  to  the  body  temperature  at  the  wall.  The  signifi¬ 
cance  of  the  condition  on  heat  flux  q^  xvi 1 1  be  considered  in 
the  next  section. 

The  preceding  development  was  based  on  laminar 
f 1 ow .  The  similarity  requirements  are  assumed  to  be  the 
same  for  turbulent  flow.  In  this  connection,  the  surface 
roughness  of  the  model  and  the  turbulence  level  of  the  tunnel 
should  be  kept  small  to  conform  with  the  prototype,  and  to 
help  locate  transition  at  a  somewhat  similar  location. 


2 . 2  Heat  _  Conduct  ion  S  irni  lar  i ty 

The  primary  means  of  heat  flow  in  the  interior  of 
the  body  is  by  conduction.  The  conduction  of  heat  is  char¬ 
acterized  by  the  classical  Fourier  equation 


**  (  JT ) 


<L  f  K  2  T  )  =.  p 


C  2L 

P  ’ZjT 

9-t 


(2.15) 


This  condition  is  automatically  fulfilled  by  a  scaled  model. 
However,  it  also  requires  that  any  subsequent  structural 
deflections  also  be  in  scale. 
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where  K,  Cp,  and  pg  are  the  heat  conductivity,  specific  heat, 
and  density,  respectively,  of  the  body  material.  Introducing 
the  non-dimensional  parameters 


the  heat  conduction  equation  becomes 


(2.16) 


(2.17; 


Here,  the  resulting  non-dimensional  parameters  that  must  be 
equal  for  both  model  and  prototype  are 

!C,U  K  ,  Cr 

Ll  (2.18) 

where  *»  Ko/PBCpc  is  the  thermal  diffussivity  of  the  material. 

The  first  parameter  is  the  Fourier  number  and  serves  to  define 

the  characteristic  time  of  the  heat  conduction  process.  The 

«=  ** 

latter  two  parameters,  K  and  C.5,  are  automatically  simulated 
if  the  same  materials  and  same  temperatures  arc  utilised  for 
both  model  and  prototype.  If  this  is  not  the  case,  some  re¬ 
laxation  in  true  similarity  will  generally  result,  although 
this  may  not  be  very  important  in  certain  cases  (quite  often, 
the  assumption  of  a  constant  K  and  Cn  yields  good  results  for 
not  too  large  temperature  ranges). 

The  boundary  condition  for  the  heat  conduction  prob¬ 
lem  is 
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(2.19 


at  the  body  surface  (see  also  Eq.  (2.14)).  In  the  above,  the 
right  hand  side  represents,  respectively,  the  heat  input  from 
the  aerodynamic  heating,  and  the  heat  loss  due  to  radiation 
from  the  body.  Non- dimensional! zing  this  as  before  yields 


Kdl)  ~(K\(k  if)  ,  (€•  To  i  )  er  T 

^ n  h»r  i  K»J(  ^ hn  J 


(2. 20) 


The  non-dimensional  parameters  that  must  be  equal  for  both 
model  and  prototype  for  simulation  of  the  boundary  condition 
are 


km  t  To  L  A,  K ,  e 


(2.21) 


The  non-dimensional  temperature  gradient  ( <2  ¥/ £  n)  °f  Eq. 
(2.20)  must  be  the  same  for  both  model  and  prototype.  This 
will  be  assured  if  the  previous  similarity  conditions  on  the 
aerodynamics  are  fulfilled. 

The  first  condition  of  Eq.  (2,21)  expresses  the 
ratio  of  heat  flux  from  the  boundary  layer  to  the  heat  con¬ 
ducted  into  the  body  and  the  second  condition  expresses  the 
ratio  of  heat  lost  by  radiation  to  heat  conducted  into  the 
body.  The  k<»/Ko,  and  K  conditions  will  be  automatically 
satisfied  if  the  same  materials,  gases,  and  temperatures  are 
utilized  for  the  model  and  prototype.  If  not,  some  relaxation 
of  the  similarity  conditions  as  discussed  before  will  result. 
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Other  boundary  conditions  such  as  at  an  insulated 
wall  (K  ~  0)  will  not  introduce  any  additional  similarity 
parameters , 

For  heat  conduction  similarity  then,  the  model  and 
prototype  must  be  geometrically  similar,  the  temperature  dis¬ 
tribution  Tair  must  be  similar,  and  the  parameters  of  Eqs. 

(2.13)  and  (2.21)  must  be  the  same  for  both.  If  radiation  is 
assumed  negligible,  and  if  model  and  prototype  are  of  the  same 
materials  and  temperatures  and  in  the  same  gas  mediums,  the 
heat  conduction  similarity  requirements  reduce  simply  to 
maintaining  the  same  Fourier  number  tC„t0/L^.  Thus  at  corre¬ 
sponding  non-dimensional  times  the  non-dimensional  temperature 
T  will  be  the  same  for  model  and  prototype  at  corresponding 
points  in  the  body. 

2 .3  St re ss  and  Deflecti on  Similarity 

The  equations  governing  the  stresses  and  deflections 
of  a  three-dimensional  heated  elastic  body  can  be  written  in 
terms  of  three  equilibrium  and  six  stress-displacement  equa¬ 
tions.  The  first  of  the  three  equilibrium  equations  is 

d  (fj  t  d'f+i  +  JTy?  =.  P  9  v*  -  p  (2.22) 

dl  d*  d]  3  9^ 


where  the  terms  on  the  right  hand  side  represent  the  inertial 
loading  and  gravity  loading,  respectively.  The  first  of  the 
six  stress-displacement  equations,  allowing  for  large  struc¬ 
tural  deflections,  is 


jC 

c  V 


du 

e?y 


[Yd*/  W/  i^  /  |  £  [ 


*  AT 


(2.23) 
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All  the  remaining  equations  are  similar  in  form  to 
equations  (2.22)  and  (2.23).  In  the  above,  Pg  is  the  density 
of  the  material,  E  is  the  modulus  of  elasticity,  y  is  Poisson's 
ratio,  e<  is  the  coefficient  of  thermal  expansion,  and  A  T  is 
the  temperature  rise  above  some  fixed  initial  value.  The  un¬ 
knowns  are  the  six  stresses  <T  x,  c Ty,  V  z,  7xy,  X  xz,  Tyz, 
and  the  three  displacements  u,  v,  v/.  The  above  equations  are 
nonlinear  in  the  displacements  u,  v,  w,  but  it  has  been  as¬ 
sumed  that  the  stress-strain  relationship  is  linear. 

If  now  the  following  non-dimensional  parameters  are 

introduced 


<T*  = 

<r> \ 

77  ' 

^ \ 
il 

<r; 

T  =»  -X- 

T. 

Txy  as 

c n  ' 

-  Ifi. 

j  =»  3k/. 

n 

7  =  -6 
'  77 

7  = 

_y 

a  / 

,  ?  -  , 

_ 

46  as.  - 

a* 

(2.24) 

AT  ~ 

/ 

Uc 

iv  =•  j£i  , 

e  =  £-,  5 

.  s 

SJ  ^  3  V 

'  i7 

the  above  equations  become 


Equations  (2.23)  and  (2.26)  introduce  the  following  non-dimen¬ 
sional  parameters 

Pi  U‘L  t  kL.  t  U.T.L 

C77  io*  £.  '  U» 

£  ,  5  ,  V  , 

The  non-dimensional  temperature  distribution  T  of  Eq.  (2.26) 
must  be  the  same  for  model  and  prototype.  This  will  be  assured 
if  the  previous  similarity  conditions  on  the  aerodynamics  and 
heat  conduction  are  fulfilled. 

The  first  condition  of  Eq.  (2.27)  represents  a  con¬ 
dition  on  dynamics  (for  static  problems  this  is  neglected), 
while  the  second  and  third  conditions  represent  the  major 
requirements  on  static  stress-deflection  similarity.  The 
parameter  a*/L  requires  the  displacements  of  the  structure  to 
be  in  the  same  scale  as  the  characteristic  dimensions  of  the 
body  and  is  important  when  considering  large  deflections.  If 
u*/l  is  maintained  separately  then  the  first  three  parameters 
of  Eq.  (2.27)  are  reduced  to 


&  , 


<r. 


(2.27) 


/ 


y 


9  To 


(2.23) 


However,  if  the  assumption  of  small  deflections  can  be  made, 
the  second  term  of  the  left  hand  side  of  Eq.  (2.26)  may  be 
neglected  and  the  parameter  u* /L  no  longer  appears  by  itself. 
Thus,  under  small  deflections,  the  structural  deformations 
need  not  be  scaled  for  stress-deflection  similarity,  and  the 
parameters  of  Eq.  (2.28)  will  appear  in  combination  with  w* II, 
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as  in  Eq.  (2.27).  The  fifth  parameter  of  Eq.  (2.27)  arises 
from  the  gravitational  forces.  The  last  four  conditions  of 
Eq.  (2.27)  are  automatically  satisfied  if  the  came  materials 
at  the  same  temperature  are  used.  Otherwise,  as  before, 
some  relaxation  of  similarity  will  result. 

For  this  elastic  body,  the  boundary  conditions  at 
the  surface  can  be  stated  either  in  terms  of  prescribed  dis¬ 
placements  or  prescribed  normal  and  tangential  stresses. 

Either 

it  -  (u  )  prescribed,  /%?-■=■  ^i/^prescribed,  vV  --  (w)  prescribed 
or 


<Tn  -  )  prescribed  -  a  'fpr 

\Jt  -  ^prescribed 


(2.29) 


The  displacement  boundary  conditions  do  not  intro¬ 
duce  any  new  parameters.  The  normal  stress  boundary  condi¬ 
tion  above  (in  which  pA  represents  the  Imposed  aerodynamic 
pressure  distribution,  and  pp  any  additional  applied  non- 
aerodynamic  loading),  is  generally  the  principal  boundary 
condition  and  it  can  be  non-dimensionali^cd  to  give. 


rn 


'ft*  Fa  + 


(2. 30) 


where  pA  °  PA/ .  Similarity  requires  that  the  parameter 
(T  ^  given  above  must  be  equal  on  both  model  and  prototype. 
I  he  p^  term  appearing  in  cr"n  will  itself  be  similar  if  the 
previous  similarity  conditions  on  the  aerodynamics  are  ful¬ 
filled. 


Mote  however,  that  consideration  of  aerodynamic  similarity 
(see  section  2.1),  requires  generally  that  W*/L  be  scaled. 
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It  is  to  be  noted  that  the  similarity  condition  on 
<T  n  involves  the  sum  of  pA  (p^  /(To)  and  (pF/do) .  Since  gener- 
ally  the  aerodynamic  pressure  distribution  ”pA  is  not  known  in 
advance,  particularly  in  aeroelastic  type  models,  the  mainte¬ 
nance  of  ?n  similarity  requires  that  the  parameters 


(fc 


(2.31) 


must  each  be  maintained  similar  individually  (the  already 
being  maintained  from  aerodynamic  similarity) .  However,  in 
certain  non-aeroelastic  type  tests  where  pA  can  possibly  be 
estimated,  similarity  can  be  achieved  through  Eq.  (2.30) 
without  the  more  restrictive  requirements  of  Eq.  (2.31). 


Assuming  that  the  temperature  distribution  T  in  the 
body  and  the  aerodynamic  pressure  distribution  pA  have  been 
simulated,  the  parameters  for  stress  and  deflection  similarity 
are  given  by  Eqs.  (2.27)  and  (2.30). 


2.4  Summary  of  General  Similarity  Requirements 

The  general  similarity  parameters  for  aerothermo- 
elastic  modeling  derived  in  the  proceeding  sections  may  be 
summarized  as  follows 

a)  External  Aerodynamic  Flow 


Up 

>  L 


(2 .  ~>2) 


b)  Heat  Conduction 


,  Koto  t  6o  ToL  tK,CF/  e,  A 

Ko  TF~  K. 


(2.33) 
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c)  Stress  and  Deflection 


<r0L 

aio  L>  Uo  ^ 

7%  +- 

£e  Ho 

_____  2T  '  — 

or 

Pq  Us  Li 

)  Ps  ,  EL  ,  ct; 

V  j 

Vo 

(2 

0~o  ioX 

(To 

These  parameters  are  for  the  most  general  time 
dependent  aerotherrnoelastic  problem.  They  have  been  derived 
subject  to  the  assumptions  stated  previously,  and  apply  gener¬ 
ally  to  elastic  bodies  below  about  1000°F  (i.e.,  M  less  than 
about  3.5). 

It  is  possible,  by  suitable  combinations,  to  alter 
some  of  the  above  similarity  parameters  into  more  compact 
and  familiar  forms.  Theu*/L  condition  of  Eqs.  (2.32)  and 
(2.34)  can  be  introduced  into  the  1st,  2nd  and  5th  parameters 
of  Eq.  (2.34),  to  result  in  CT0/E0,  <*0To,  and  PBL2/cr0t02 
respectively.  The  (p„  pA  +  pp ) /<T0  parameter  of  Eq.  (2.34), 
as  discussed  in  section  2.3,  implies  the  maintenance  of  the 
separate  parameters  p*  /<To  and  p F/<r0.  The  p*  /(To  parameter, 
which  then  may  serve  to  define  the  reference  stress  level 
(To,  can  be  combined  with  the  <T0/E0  to  yield  the  basic  aero- 
elastic  similarity  parameter  prf>/E0.  This  basic  aeroelastic 
parameter  expresses  the  ratio  of  aerodynamic  forces  to  stiff- 
ness  forces.  The  parameter  Pgl//<r0t0z  can  also  be  transform¬ 
ed,  by  combination  with  the  similarity  conditions  Ut0/L, 


It  could  also  have  been  expressed,  by  further  combination 
with  the  perfect  gas  relation  p^  «  P»  U2/^m2  and  the 
and  H  similarity  relations,  in  the  alternate  form  q/E0l 
where  q  is  the  dynamic  pressure. 
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WZ-'-  -  W- 


f 

f 

2  2 

P*>  /(Tq>  <T >  30(3  the  perfect  gas  relation  pw  -  P«,  U  /^M  , 

into  the  simple  form  Pg/P«»  ,  which  is  the  mass  density  ratio 

parameter  important  in  flutter  work.  Similarly.,  the  parameter 

f>*tL/  Oo  can  be  changed,  by  combination  with'  PnLZ /<TntJ'  and 

rj*  /  ,  "  ° 

Ut0A>  into  the  form  This  may  be  identified  with  the 

Froude  number,  which  expresses  the  ratio  of  gravity  forces 
to  aerodynamic  forces. 

The  general  aerothermoelastic  similarity  parameters 
of  Eqs.  (2.32),  (2.33),  and  (2.34)  may  be  rewritten  in  terms 
of  these  above  combined  alternate  forms  as. 


M,  Re,  fz.  ,  L  ,  K  ,  *.T.,  Pr 

£-  Ko 

y ,  ~  ~  t  t  ^  ^  / *  >  m  >  ^  /  Cp 

Ko  fL 

i,  *  ,  v ,  e,  u&l,  Koto  t  ,  ZL 

^  L7'  L*  (To  (To  7^ 


(2,35) 


The  first  six  parameters  of  the  above  are  the  principal  con¬ 
ditions  to  be  satisfied  for  aerothermoelastic  similarity. 

The  next  five  are  generally  somewhat  weaker  requirements. 

The  next  nine  double-barred  quantities  reflect  the  variations 
in  the  parameters  for  conditions  departing  from  the  reference 
temperature  T0.  For  same  model  and  prototype  materials  at 
the  same  temperatures,  these  are  all  identically  satisfied. 
The  remaining  six  parameters  serve  to  define  the  reference 
values  o£  to  ,  Uo  ,  (To  and  7^  to  be  used  in  the  non-dimen- 
sionalizations.  For  the  case  of  static  aerothermoelastic 
models,  the  parameters  F*/Peo  s  ^°/u  can  be  dropped  in  the 
above,  and  the  Froude  number  parameter  U fyb  will  be  replaced 
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by  the  original  fa  ty<To  parameter  of  Eq.  (2.34)  . 

To  investigate  the  application  of  the  similarity 
parameters  of  Eq.  (2.35),  consider  an  experiment  in  which  a 
scaled  model  is  constructed  of  the  same  material  and  tested 
in  the  same  gas  and  at  the  same  free  stream  temperature  as 
the  prototype.  It  is  assumed  that  the  model  will  represent 
a  static  aerothermoeiastic  test  and  that  the  effects  of  radi¬ 
ation,  gravity  forces,  and  all  non- aerodynamic  loads  fF  on 
both  model  and  prototype  are  negligible.  The  model  scale 
factor  is  designated  by 


h 


(2.36) 


where  the  subscript  P  denotes  a  prototype  quantity  and  M  a 
model  quantity. 

Under  the  above  stated  conditions,  the  similarity 
parameters  of  Eq.  (2.35)  reduce  to. 


M,  Re  ,  fn,  ,  Ko  tc  , 

L'  ~L  IT 

The  following  relations  will  be  obtained  between  model  and 
prototype  quantities  by  maintaining  all  the  above  similarity 
conditions,  exceot  for  the  fm/se  condition. 
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(2.38) 


7f  (  Uo)r 


ihe  remaining  'f«> /Eo  condition,  which  arises  from  the  aero¬ 
dynamic  forces  on  the  body,  will  result  in  the  relation, 

(f'O*  —  (fm)p  (2 

This  condition  contradicts  the  fifth  condition  of  Eqs.  (2.38) 
which  arose  from  consideration  of  the  Reynolds  number  condi¬ 
tion.  Hence  for  this  simplified  case  of  similar  materials, 
similar  gas,  and  similar  temperatures  in  a  static  aerothermo- 
elastic  test,  conflicts  arise  among  the  similarity  parameters, 
and  one  is  unable  to  correctly  simulate  the  prototype  situa¬ 
tion.  This  is  also  pointed  out  in  Refs.  3,  5,  6,  7,  8  and  10. 

Three  other  conflicting  requirements,  though  gener¬ 
ally  less  severe,  would  also  have  arisen  had  the  more  general 
dynamic  aerothermoelast ic  case  been  considered.  The  first 
would  have  been  a  conflict  in  the  characteristic  times  t0 
resulting  from  the  U*m/b  and  conditions  for  same 

materials.  The  former  would  have  given  a  ft*)* /$*),  ratio 
of  1/n  rather  than  the  1/n2  of  the  latter.  This  is  not, how- 
ever,  a  serious  conflict  because  of  the  substantially  differ¬ 
ent  time  constants  associated  with  the  vibrational  and  aero¬ 
dynamic  phenomena  (very  short)  and  those  of  the  thermal  heat¬ 
hs  phenomena  (quite  long) .  Only  if  these  are  of  about  the 
same  order  of  magnitude  and  if  mutual  coupling  is  possible 


.39) 
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between  these  two  phenomena  would  it  be  necessary  to  have  a 
single  common  reference  time  to  for  the  specific  aerothermo- 
elastic  test,  A  second  conflicting  requirement  might  arise 
from  the  radiation  Ko  condition.  For  these  similar 

materials  at  the  same  temperatures,  the  scale  ratio  would 
prevent  satisfaction  of  this.  Unless  rather  high  temperatures 
are  involved  though  this  radiation  would  tend  to  be  unimpor¬ 
tant.  Even  then  it  might  be  alleviated  to  some  extent  by 
altering  the  surface  characteristics  and  hence  the  emissivity 
6  of  the  model  somewhat.  Another  source  of  conflict  lies  in 
the  Froude  number  V  / jL  .  in  most  aeroelastic  problems 
though  this  effect  is  small.  The  possible  conflict  involving 
the  mass  density  ratio  P&  /fco  can  be  alleviated  somewhat  by 
adding  weights  artificially  to  built-up  type  models. 

However,  the  fundamental  conflict  in  similarity 
remains  that  concerning  the  Mach  number  M,  Reynolds  number 
Re,  and  aeroelastic  parameter  / Eo  >  and  means  of  getting 

about  this  must  be  investigated. 

One  means  of  doing  this  is  to  attempt  to  use  differ¬ 
ent  model  and  prototype  materials,  in  different  gases,  and  at 
different  temperatures.  This  approach  is  studied  in  Section  3. 

Another  approach  is  to  develop  less  restrictive 
similarity  parameters  by  considering  more  specialized  situa¬ 
tions  than  the  very  general  ones  given  here  for  the  aerody¬ 
namics,  heat  conduction,  and  stress-deflections.  For  these 
specialized  situations  it  is  possible  to  develop  similarity 
laws  permitting  even  geometrical  distortions  of  the  scaling 
down  process.  This  approach  which  has  been  successful  in 
developing,  for  example,  cold  flutter  model  techniques,  is 
explored  in  Section  4. 

A  further  approach  to  this  conflict  is  to  relax 
one  of  the  requirements  for  similarity.  Either  the  M,  the 
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Re,  or  the  fa/Sc  requirement  can  be  left  unsatisfied.  Con¬ 
sideration  of  the  specialized  similarity  situations  mentioned 
above  may  help  give  a  better  clue  as  to  what  parameters  are 
important  in  considering  specific  types  of  phenomena,  and 
v/hat  parameters  may  therefore  be  dropped.  For  example,  the 
Re  requirement  can  be  shown  to  be  associated  mainly  with  the 
aerodynamic  heating,  and  if  this  heating  is  unimportant,  the 
Re  requirement  can  be  dropped.  This  is  effectively  what  is 
done  in  ordinary  cold  aeroelastie  testing.  Or,  if  the  aero¬ 
dynamic  loadings  and  subsequent  structural  distortions  were 
not  too  important,  the  ‘fm/£o  requirement  could  be  disregard¬ 
ed.  This  is  done  in  ordinary  heat  transfer  and  friction 
studies  on  models.  Such  types  of  "restricted  purpose"  model 
tests  can  be  of  use  in  investigating  carefully  delineated 
problems.  The  consequences  of  relaxing  some  parameters  are 
discussed  in  Section  5. 

It  is  of  interest  to  note  that  another  means  of 
resolving  this  conflict  between  M,  Re,  and  fm /So  is  possible 
in  certain  non- aeroelastie  cases  by  placing  additional  non- 
aerodynamic  surface  loads  on  the  model.  As  discussed  in 
Section  2.3,  the  normal  stress  similarity  parameter  is  more 
generally  .  Using  this  relation,  the  separate  f*/So 
parameter  of  Eq.  (2.35)  will  be  replaced  by  the  combined 
requirement ? 

fm  f*  +  (2.40) 

So 

Appling  this  to  a  prototype  without  additional  loads  and 
a  model  with  additional  loads,  gives 


(2.42) 


This  will  result  in  the  requirement  of  placing  additional 
external  loadings  (ff)^  over  the  model  of  magnitude 

where  is  the  presumed  known  aerodynamic  pressure  distri 
bution.  In  an  aeroelastic  type  model,  the  determination  of, 
and  the  supplying  of  this  additional  surface  loading  might 
present  considerable  difficulty* 
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3.  CONSIDERATION  OF  DIFFERENT1  MATERIALS  AND  DIFFERENT  GASES 
3,3  Different  Materials 

In  the  preceding  section  the  similarity  parameters 
of  the  general  aerothermoelastic  problem  were  derived  and 
some  of  the  inherent  conflicts  leading  to  the  requirement  of 
length  scale  ratios  of  1  were  indicated.  It  is  of  interest 
to  investigate  the  use  of  different  model  and  prototype  mate¬ 
rials  and  gases  other  than  air  for  the  test  medium  in  match¬ 
ing  these  parameters.  With  this  approach  model  stream  tem¬ 
peratures  need  not  also  equal  prototype  stream  temperatures 
and  length  scale  ratios  greater  than  1  are  generally  possible. 

The  parameters  required  for  similarity  are  given 
in  Eq.  (2.35).  The  principal  similarity  parameters  to  be 
matched  here  will  be  the  M,  Re,  f*  / Eo  ,  A^/Ko  ,  and  o(oTo 
conditions.  In  attempting  to  match  these,  it  is  convenient 
to  transform  the  Re  requirement  by  introduction  of  the  per¬ 
fect  gas  law,  Eq.  (2.5),  and  combining  with  the  M  and  fk/£o 
requirements  into  the  new  parameter, 

L  Eo  riK  o 

■J7J  r, 

From  the  above,  it  can  be  seen  that  mutual  satisfaction  of 
the  M,  Re,  and  -fa>/E0  conditions  will  result  only  for  scale 
ratios  of 


Thus  to  achieve  scale  ratios  greater  than  1  it  is  desireable 
for  the  model  material  to  have  a  higher  modulus  of  elasticity 
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than  the  prototype  and  to  use  a  gas  heavier  than  air  in  the 
tunnel  at  a  lower  temperature  than  the  prototype. 

In  addition  to  satisfying  M,  Re,  /e*  ,  /k<> 

and  <><o  lo  t  the  density  ratio  fa /fa  should  also  be  satisfied 
(for  dynamic  tests  only)  as_  well  as  the__non-dimensional  mate¬ 
rial  and  fluid  properties  n  ,  Cr  ,£,<*,  V  ,  It  ,  . 

The  Pg  /fa  parameter  can  often  be  manipulated  by  addition 
of  weights  to  built-up  type  structures.  In  satisfying  the 
latter  double-barred  parameters  one  is  essentially  matching 
the  slopes  of  model  and  prototype  materials  at  the  reference 
temperatures  C7®%i  and  (' 7*)p  .  This  is  extremely  difficult 

and  the  problem  may  be  facilitated  by  selecting  T0  to  lie 
somewhere  in  the  neighborhood  of  the  expected  maximum  temper- 
ature.  In  this  manner  K  ,  ,c*  etc.  are  approximately 

one,  in  this  neighborhood,  the  error  depending  on  the  differ¬ 
ences  in  slopes  of  material  properties.  Of  course,  when  the 
same  materials,  gas,  and  temperatures  are  utilized  these 
quantities  are  automatically  satisfied.  It  should  be  noted 
that  if  a  material  or  fluid  property,  ~)\  ,  can  be  represented 
by  a  temperature  dependence  of  the  form 

}  ~  T  __  (-.) 

for  some  temperature  range,  then  will  be  automatically 
satisfied  since  ^"/V  is  the  same  for  model  and  prototype  for 
temperature  similarity.  The  radiation  condition  and  Froude 
number  will  be  assumed  to  be  negligible  here. 


The  fluid  properties  A  ,E,£r,  which  are  referred  to 
free  stream  temperature  rather  than  the  reference  tempera¬ 
ture  T n,  are  not  1.  However,  these  properties  are  repre- 
sentible  in  the  form  of  Eq.  (3.3)  for  this  temperature 
range  (0  to  1000°F)  and  hence  these  quantities  are  also 
generally  satisfied. 


ASRL  TR  75-1 


25 


-  II  I  I  1 

***•*••' 


' 


■ 


Five  structural  materials,  aluminum,  magnesium, 
stainless  steel,  tncone l=x, and  titanium  have  been  considered 
wi tn  respect  to  being  used  either  as  prototype  or  model  mate¬ 
rials.  Material  properties  for  these  are  given,  for  example, 
in  Refs.  26,  27  and  28.  For  these  materials  and  for  air  as 
the  wind  tunnel  gas,  the  parameters  &d>/Ko  ,  t  and  Eq. 

(3.2)  have  been  plotted  as  a  function  of  T  in  Figs.  1,  2. 
and  3.  These  plots  illustrate  concisely  the  difficulties  en¬ 
countered  in  satisfying  the  parameters  for  general  similarity. 
Figure  2  shows  that  similar  materials  will  satisfy  the  c*eT0 
condition  only  if  model  and  prototype  temperatures  are  the 
same.  However,  Fig.  3  shows  that  if  similar  materials  and 
similar  temperatures  are  utilized  the  scale  ratio  must  be  1. 

Thus  to  simulate  the  To  condition  along  with  M,  Re  and  fa / Eo 
for  scale  ratios  greater  than  1,  model  and  prototype  materials 
and  temperatures  must  be  different.  Figure  1  shows  that  the  kmjKt 
condition  cannot  be  satisfied  if  aluminum  or  magnesium  are 
used  for  matching  with  either  inconel-x,  stainless  steel,  or 
titanium.  With  a  small  degree  of  error  inconel-x,  stainless 
steel  and  titanium  may  be  used  for  matching  with  each  other. 

Some  examples  of  the  possible  scale  ratios  resulting 
from  the  use  of  different  materials  and  utilizing  air  as  a 
tunnel  gas  are  given  in  Fig.  4.  It  is  apparent  that  mutual 
satisfaction  of  M,  Re,  fa / *£•  ,  k»/Ko  and  <*•  7S  is  at  best 
only  approximate  for  certain  particular  model-prototype  com¬ 
binations  and  even  then  for  scale  ratios  not  very  far  differ¬ 
ent  from  1.  For  example,  the  best  combination  appears  to  be 
a  titanium  prototype  with  a  stainless  steel  model  at  a  scale 


koo  and  p, <*>  have  been  replaced  here  by  k0 
mentioned  above,  both  these  can  generally 
the  form  of  Eq.  (3.3)  in  this  temperature 


and  p,  since  as 
be  expressed  in 
range . 
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ratio  of  3. 
condition, 
less  steel. 


4.  An  aluminum  prototype,  because  of  the  !k0 
cannot  be  simulated  by  a  model  constructed  of  stain 
incopel-x  or  titanium  (see  Fig.  3). 


It  should  be  noted  that  similitude  of  thermal  stresses 
and  deformations  depends  on  both  the  temperature  distribution 
T  and  the  parameter  <Ae~fo  as  mentioned  in  section  2,.  Thus, 
since  simulation  of  T  depends  on  *m/Ko  ,  simulation  of  thermal 
stresses  will  depend  on  as  well  as  on  .  Further¬ 

more,  since  the  temperature  ratio  T  is  the  same  for  model  and 
prototype,  the  initial  temperature  of  the  model  must  be  in  the 
same  ratio  to  the  prototype  initial  temperature  as  (~T°)  M  /(*;,  • 
Thus  if  is  less  than  (7 °)p  it  may  be  necessary  to  refri- 

geruie  the  model  initially.  If  the  model  is  of  uniform  mate¬ 
rial.  througnout  this  presents  no  problem.  However,  if  various 
portions  of  the  model  are  constructed  of  different  materials 
then  stresses  may  arise  as  a  result  of  unequal  coefficients 
of  expansion  of  the  different  materials. 


3.2  Use  of  Gases  Other  Than  Air  in  the  Wind  Tunnel 

The  use  of  a  gas  other  than  air  as  the  test  medium 
provides  another  approach  to  the  problem  of  simulating  the 
parameters  of  the  general  aerothermoelastic  problem. 

Helium  has  received  attention  as  a  substitute  for 
air  primarily  because  of  its  low  condensation  temperature, 
eliminating  the  necessity  of  heating  the  wind  tunnel  gas  (Ref. 
11) .  Helium  is  also  attractive  because  of  its  high  thermal 
conductivity.  A  model  constructed  of  either  aluminum  or  mag¬ 
nesium  and  tested  in  helium  will  match  the  condition 

with  prototypes  constructed  of  inconel-x,  stainless  steel,  or 
titanium,  with  much  less  error  than  if  air  were  used  (Fig.  1) ♦ 
However,  if  the  Mach  number,  Reynolds  number  and  the  /’<’/& 
conditions  are  also  to  be  satisfied  simultaneously  (Eq.  (3.2)), 
a  lower  scale  ratio  ^/l ^  results  with  helium  than  with  air. 
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This  is  due  primarily  to  the  lower  molecular  weight  of  helium. 
It  may  be  seen  from  Eq.  (3.2)  that  a  model,  made  of  the  same 
material  as  the  prototype,  and  tested  in  a  substitute  gas, 
will  have  a  scale  ratio  ?fL-M  equal  to  jX^tJ  time.j 

the  scale  ratio  of  the  equivalent  test  in  air  at  the  same  Mach 
number,  Reynolds  number,  temperature,  and  pressure.  Thus  the 
desirable  characteristics  of  a  substitute  gas  are  greater  mole¬ 
cular  weight  and  lower  viscosity  than  air.'  Fig.  5,  derived 
from  the  data  of  Refs.  12,  13  and  14,  presents  some  properties 
of  gases  which  have  en  considered  as  possible  substitutes 
for  air.  Regarding  scale  ratio,  helium  is  the  least  attrac¬ 
tive,  providing  a  scale  ratio  .33  times  that  available  with 
air,  while  F10,  the  heaviest  gas,  provides  a  scale  ratio 
3.86  times  greater  than  is  possible  with  air  under  the  same 
tunnel  conditions. 

The  substitute  gases  of  Fig.  5  do  not  have  the  same 
specific  heat  ratio  y'as  air.  At  the  higher  supersonic  speeds, 
the  normal  pressure  distribution  on  the  lifting  surface  is 
dependent  on  "^and  failure  to  simulate  this  parameter  is  one 
of  the  primary  disadvantages  of  using  a  substitute  gas.  Thus 
Cz  f10  with  1.06,  which  was  the  most  attractive  from  scale 
ratio  considerations,  now  appears  undesirable.  Reference  12 
shows  that  a  means  of  simulating  ^  exactly  is  possible  by 
properly  mixing  a  heavy  monatomic  gas,  such  as  argon,  krypton 
or  xenon,  with  a  suitable  heavy  polyatomic  gas.  Three  illu¬ 
strations  of  this  are  also  given  in  Fig.  5,  in  which  CBr  fj,, 

SF^  and  CF,  are  each  combined  with  argon  in  proportions  that 
yield  a  ^  of  1.4  for  the  mixture.  The  improvement  in  scale 
ratio,  for  each  of  these  mixtures,  is  seen  to  be  approximately 
1.5  times  greater  than  is  possible  for  the  corresponding  test 
in  air. 

There  are  other  advantages  to  be  realized  by  using 
a  heavy  gas  for  aeroelastic  testing.  At  a  fixed  temperature, 
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the  speed  of  sound  is  lower  in  a  heavier  gas,  consequently, 
tests  conducted  at  a  given  Mach  number  and  temperature  in  a 
heavy  gas  will  be  run  at  a  lower  velocity  than  the  correspond¬ 
ing  test  in  air.  The  advantage  of  a  reduced  velocity  enables 
the  Froude  number  UVfL  to  be  satisfied,  if  this  is  important, 
for  scale  ratios  other  than  1.  Also,  oscillation  frequencies 
in  a  flutter  test  will  be  reduced,  and  more  important,  the 
horsepower  required  of  the  tunnel  is  only  a  fraction  of  that 
required  for  air  at  the  same  tunnel  conditions,  (see  Ref.  12). 

By  using  a  combination  of  both  different  gases  and 
different  materials  at  different  temperatures,  somewhat  larger 
scale  ratios  can  be  achieved.  Combining  the  results  above 
with  those  of  section  3.1,  scale  ratios  of  approximately 
3.4  x  1.5  ~  5  can  be  roughly  achieved. 

It  has  been  shown  in  this  section  that  by  using 
different  materials  and  gases  satisfaction  of  the  similarity 
parameters  of  the  general  aerothermoelastic  problem  is  ex¬ 
tremely  difficult.  Where  scale  ratios  other  than  1  are  pos¬ 
sible  they  are  still  not  very  large,  and  appear  inadequate 
for  testing  complete  wings  or  wing-body  combinations. 
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4.  SPECIALIZATIONS  OF  ^EROIHERHOELASTIC  SDlILARITY,ijWS 

The  scaling  requirements  of  the  general  aerothermo- 
elastic  problem  considered  in  the  previous  sections  have  been 
shown  to  be  very  stringent  and  difficult  to  satisfy  completely* 
It  is  possible  however.,  that  by  considering  specific  aerother- 
moelastic  problems  and  hence  using  more  specialized  governing 
equations  than  the  general  ones  of  the  previous  sections,  to 
arrive  at  less  restrictive  similarity  conditions.  This  tech¬ 
nique  may  also  help  point  out  the  significance  of  the  various 
general  similarity  conditions. 

4,1  Specialized  Aerodynamic  Similarity 

The  boundary  layer  concept  introduced  by  Prandtl 
has  permitted  the  solution  of  aerodynamic  flow  about  an  object 
to  be  divided  into  two  separate  problems;  the  first  dealing 
with  the  de' ermination  of  the  pressure  distribution,  and  the 
second  dealing  with  the  determination  of  the  viscous  and  heat¬ 
ing  effects  in  a  thin  boundary  layer  region  close  to  the  body. 
This  separation  appears  valid  for  high  Reynolds  numbers  and 
for  M  less  than  about  5,  and  is  widely  used  in  many  circum¬ 
stances  for  estimating  pressure  distribution  and  neat  transfer 

data. 

a)  Aerodynamic  Pressure  Distribution 

For  determining  the  pressure  distribution  about  an 
object,  the  governing  equations  are  assumed  those  of  a  non- 
viscous,  non-heat -conducting  fluid.  Under  these  conditions, 
the  viscosity  (j,  and  heat  conduction  k  are  assumed  zero  in  the 
general  aerodynamic  equations  (Eqs.  (2.1)  to  (2.5)),  and  the 
required  similarity  reduces  to  constant  values  of 


Ut/ 

The  term,  /u  ,  arises,  as  mentioned  previously, 
from  the  boundary  condition  of  no  flow  through  the  body  and 
here  requires  the  ratio  of  the  thickness  (and/or  displacement) 
Uo  to  the  scale  length  L  ,  to  be  the  same  on  both  model  and 
prototype.  The  similarity  conditions  expressed  above  are  far 
less  stringent  than  those  of  Eqs.  (2.13)  chiefly  because  of 
the  absence  of  the  Reynolds  number  Re  requirement,  and  indeed 
are  in  wide  use  in  determining  aerodynamic  pressures,  stabil¬ 
ity,  and  flutter  data  from  wind  tunnel  models. 


In  certain  circumstances,  the  similarity  conditions 
of  Eq.  (4.1)  can  be  reduce  1  still  further.  For  small  pertur¬ 
bations  in  a  shockless  flow  at  Mach  numbers  greater  than  about 
2.5,  aerodynamic  piston  theory  (see  Ref.  15)  can  frequently 
be  used  to  describe  the  pressure  differential  between  upper 
and  lower  surfaces  of  a  thin  lifting  surface.  This  can  be 
written  in  the  form. 


—  .  "J 

/  j-Q'+O  M<dJL 

dw 

+  /  d  w 

j 

d* 

U  0-t  _ 

(*  -2) 


where  q  represents  the  dynamic  pressure,  S  represents  the 
thickness  of  the  surface,  and  w  the  displacement  of  the  mean 
line  of  the  surface.  Non-dimensionalizing  as  before,  the 
above  reduces  to. 


A  f  = 


■(MU.) 

r*  231 

/  +  ft' 

'  d*_  +(L_\dw_~ 

'  Ls  J 

.  *  1  j*  _ 

\utt)  di. 

(4.3) 


The  similarity  requirements  are  then  simply 


Af  M  $  o  *  y  f  Uto  }  S 

^  L.  ~~C~ 


(4.4) 
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It  can  be  seen  here  that  Mach  number  similarity  is  no  longer 
required  by  itself,  but  rather  in  the  combined  parameters 
M  g*  and  &  .  If  the  thickness  distribution,  3  s  */So, 
and  downward  deflection  distribution  WA.  are  similar,  the 
similarity  requirements  above  can  permit  different  Mach  num¬ 
bers  in  model  and  prototype,  provided  only  that  both  model 
and  prototype  are  adequately  characterized  by  piston  theory. 

Aerodynamic  similarity  relations  f.ox  the  same  family 
of  profiles  and  shapes  can  likewise  be  derived  in  the  subsonic, 
transonic,  and  low  supersonic  regimes.  These  constitute  the 
well  known  aerodynamic  similarity  rules  (see  for  example  Ref. 

16)  and  are  useful  in  relating  the  pressure  characteristics 
of  similar  families  of  slender  bodies. 

*  . 

In  considering  hypersonic  Mach  numbers  (M>5), 
there  is  an  interaction  between  thickness  of  the  boundary 
layer  and  the  external  flow  field  which  affects  the  external 
pressure  distribution.  This  specific  similarity  problem  is 
discussed  in  Ref.  17,  and  it  is  shown  there  that  an  addition¬ 
al  combined  parameter  M^ / J Re' should  be  included  in  the  aero¬ 
dynamic  pressure  similarity  conditions  of  Eqs.  (4.4)  together 
with  Pr,  T* /  f0  ,  and  Zp  requirements. 

In  the  vicinity  of  noses  of  blunt  bodies  at  M  “2.5, 
the  Mach  number  disappears  altogether  from  the  required  simi¬ 
larity  parameters,  as  then  Newtonian  type  flow  considerations 
take  over  (see  Ref.  18).  The  pressure  distribution  can  then 
be  expressed  simply  as 

_£  -  COS  9  (4.5) 

f'"b 


*  This  is  somewhat  greater  than  the  regions  here  investigated 


where  &  represents  the  flow  deviation  angle  at  the  surface 
of  the  body.  This  is  one  of  the  features  that  makes  shock 
tube  investigations  at  high  Mach  numbers  useful. 


b)  Viscous  and  Heating  Effects 

For  obtaining  viscous  and  heating  effects,  one  must 
look  at  the  complete  equations  as  given  by  Eqs.  (2.1)  to  (2.5). 
These  can  be  simplified  further  through  Prandtl's  order  of 
magnitude  argument,  to  yield  solutions  in  the  boundary  layer 
region  close  to  the  body.  The  required  similarity  parameters 
however  still  remain  the  same  as  given  by  Eq.  (2.13)  for  the 
most  general  conditions. 


Under  certain  circumstances,  these  parameters  can 
be  made  less  restrictive.  The  unsteady  effects  associated 
with  heating  are  generally  small,  and  hence  the  parameter 
Vto/u  can  usually  be  neglected.  For  the  two-dimensional, 
steady,  laminar  flow  over  an  isothermal,  semi- inf inite  flat 
plate,  the  boundary  layer  equations  can  be  solved  to  give  the 
following  expression  for  the  heat  flux  input 


1-  ■  ('%'L 


(4.6; 


where  ®  adiabatic  wall  temperature. 

For  thin  surfaces  with  samll  pressure  gradients  the  above 
expression  can  give  a  reasonable  estimate  to  the  amount  of 
heating  experienced  by  the  wing  in  laminar  flow. 

Upon  non-dimenslonalization  of  Eq.  (4.6)  one  obtains 
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The  above  in  combination  with  the  general  heat  conduction 
equation  at  the  surface,  Eq.  (2.19),  results  (upon  neglecting 
radiation  for  the  present)  in  the  parameters 


The  second  parameter  defines  the  reference  temperature  in 
terms  of  the  adiabatic  wall  temperature,  while  the  first  com¬ 
bined  parameter  involving  Reynolds  number  does  away  with  a 
good  many  of  the  previous  separate  aerodynamic  and  thermal 
requirements.  Thus,  the  necessity  of  similar  Mach  number  is 
not  required  for  viscous  heating  under  these  circumstances. 
(Notice,  though,  that  the  coefficient  .332  may  change  somewhat 
as  a  function  of  Mach  number,  but  this  can  readily  be  account¬ 
ed  for  in  the  above  combined  parameter).  The  presence  of 
temperature  gradients  on  the  plate  will  not  affect  the  above 
similarity  relation  if  one  assumes  the  validity  of  the  Light- 
hill  approximation  to  the  heat  transfer  (see  Ref.  19). 

For  turbulent  flows  over  thin  surfaces  with  small 
pressure  gradients,  the  analogous  similarity  parameter  to  the 
first  one  of  Eq.  (4.8)  is, 

.OZ9&  JL_  (Re)  (Pr)  (4.9) 

K0 

where  now  the  coefficient  .0296  is  somewhat  more  affected  by 
Mach  number  than  in  the  laminar  case.  As  in  the  laminar  case, 
the  presence  of  temperature  gradients  will  have  very  little 
effect  on  this  similarity  parameter  (Ref,  20). 

The  validity  of  the  above  similarity  parameters  can 
be  shown  to  extend  also  to  the  aerodynamic  heating  at  the 


"  -.i'  * 


b3W)WM 


noses  of  blunt  bodies  providing  the  constants  are  altered  to 
.763  and  .042  respectively  (Ref.  21).  The  location  of  the 
transition  point  from  laminar  to  turbulent  flow  is  however  a 
nebulous  quantity.  For  certain  circumstances  however,  it  may 
be  assumed  to  be  independent  of  Mach  number  to  a  first  ap¬ 
proximation. 

To  summarize  then,  it  is  seen  that  by  means  of  the 
Prandtl  boundary  layer  concept  for  high  Reynolds  numbers  and 
not  too  high  Mach  numbers,  that  it  is  possible  to  associate 
the  pressure  distribution  similarity  mainly  with  the  Mach 
number  parameter  M,  and  the  viscous  and  heating  similarity 
mainly  with  the  Reynolds  number  parameter  Re.  In  certain 
circumstances  (for  example,  flow  over  thin  surfaces),  it  is 
possible  to  relax  somewhat  the  requirement  of  Mach  number 
and  Reynolds  number  being  individually  similar,  and  to  obtain 
rather,  combined  parameters  such  as  Eqs.  (4.4)  and  (4.8) 
which  can  be  easier  to  satisfy.  Of  course,  if  in  any  model 
test  either  pressure  distribution  or  viscous  drag  and  heating 
is  unimportant,  the  corresponding  similarity  parameters  may 
be  neglected  without  great  loss  of  accuracy  for  the  particular 
test . 

4.2  Heat  Conduction 

a)  Thin  Plates 

In  the  heating  of  a  thin  plate,  the  temperature 
tends  to  rise  uniformly  across  the  plate  thickness  resulting 
in  a  zero  temperature  gradient  in  the  thickness  direction. 
Under  these  conditions  the  heat  conduction  equation  for  a 
thin  plate  of  thickness  <f  may  be  written  as 


where  the  last  two  terms  on  the  right  hand  side  represent 
the  aerodynamic  heating  input  and  radiation  heat  loss*  re¬ 
spectively*  over  the  surface  of  the  plate.  The  boundary 
conditions  in  this  case  need  only  be  applied  to  the  edges  of 
the  plate  and  not  to  the  surface.  Tha  conditions  for  appli¬ 
cability  of  this  thin  plate  assumption  are  generally  that 
the  quantity  AS/k  be  less  than  about  1,  where  h,  £  *  K  are 
the  heat  transfer  coefficient*  the  thickness*  and  the  thermal 
conductivity  respectively. 

Non-dimensionalizing  by  introducing  the  same  para¬ 
meters  as  in  Section  2  plus  the  additional  thickness  distri¬ 
bution  parameter  S-  8 / $0  results  in 


The  resulting  similarity  parameters  for  the  temperature  dis¬ 
tribution  in  the  thin  plate  are  then 

,  KM  ,  €.Tq3L  t  K ,  Cr ,  A ,  ?  (L 

Ko  So  L 

The  similarity  conditions  of  Eq.  (4.12)  are  less  restrictive 
than  those  of  Eq.  (2.33)  since,  by  using  a  different  material 
or  gas  for  the  model,  the  plate  thickness  So  need  not  be 
scaled  in  the  same  manner  as  the  characteristic  length  L  . 

This  is  manifested  by  the  parameter  l/k»  $*  which  replaces 
the  parameter  k»/Ko  of  Eq.  (2.21).  The  edge  boundary  con¬ 
ditions  will  introduce  no  additional  parameters  to  Eq.  (4.12). 
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If  the  further  assumption  of  no  heat  conduction  in 
the  x  and  y  directions  can  be  made,  the  left  hand  side  of 
Eq.  (4.11)  is  set  equal  to  zero  and  the  non-dimensional  para¬ 
meters  became 


io  j£a> 


3 

71  L 


A 


(<•.13) 


These  parameters  allow  even  greater  flexibility  because  there 
is  no  scaling  restriction  on  the  plate  thickness  other  than 
the  requirement  that  the  model  also  be  thermally  thin.  The 
first  parameter  of  Eq.  (4.13)  will  serve  only  to  define  the 
reference  time  tQ  of  the  thermal  response  of  the  plate. 

b)  Thick  Plates 

The  similarity  parameters  for  the  temperature  dis¬ 
tribution  in  a  thick  plate  are  the  same  as  for  the  arbitrary 
body,  Eqs.  (2.18)  and  (2.21).  However,  flexibility  in  the 
scaling  of  the  plate  thickness  is  possible  for  situations 
where  the  temperature  gradients  in  the  x  and  y  directions  are 
negligible.  Non-dimensionalizing  the  z  direction  with  respect 
to  So  instead  of  L,  and  neglecting  the  temperature  gradients 

77  ’  >  E98-  (2-15)  and  (2.19)  will  yield  the  para¬ 

meters, 

K'  L-  ,  ,  ttJLL  J  I ,  €  (4.1 

0°  do 


Therefore,  as  with  the  thermally  thin  plate  problem,  the 
reference  thickness  !•  of  a  thick  plate  need  not  be  scaled 
in  the  same  manner  as  L  if  a  different  material  or  gas  is 
used  for  the  model. 
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c)  Built-Up  Wina  Structures 

Consider  the  built-up  type  of  wing  of  Fig.  6  con¬ 
sisting  of  several  webs  and  ribs.  If  tne  skin  Is  thermally 
thin,  the  governing  non-dimensional  heat  conduction  equation 
will  be  Eq.  (4.11).  For  each  web,  the  non-dimensional  heat 
conduction  equation  is 


and  the  boundary  conditions  at  each  skin-web  junction  are 

7=  .  =F.  (4.16a) 


'w 


\  }$,*$•! 


Kw  3  ^ 


(4.16b) 


where 


JW 


Z-  ,  K«  =  K*.  ,  Cp  = 

U  Ki»  Zp 


(4 . 16c) 


y^e> 


The  resulting  non-dimensional  parameters  from  Eqs.  (4.11), 
(4.15)  and  (4.16),  are 


r  * 

L  f  t  dw0^e_  t  6;os7oL^  / 


a- 


cT<S0  ffs*  ^ 


<0 


(4.17) 


The  first  two  parameters  relate  to  the  skin  similarity,  while 
the  third  relates  to  the  web  time  history.  The  fourth  para¬ 
meter  relates  to  the  junction  boundary  condition.  The  para¬ 
meter  S /fa  requires  that  model  and  prototype  have  the  same 
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Internal  heat  flOW  PathS-  ^ite  in  regions  away  from 

he  leading  and  trailing  edges  of  the  wing,  conduction  along 
the  outer  skin  is  negligible  and  the  terms  containing  temper- 
ature  gradients  in  the  x  and  y  directions  in  Eqs.  (4.n), 

(4  15)  and  (4.16)  may  be  neglected.  Therefore,  the  parameters 

Eq.  (4.13)  will  apply  for  the  skin  and  the  similarity  para¬ 
meters  for  the  skin-web  combination  become 


Aoo 


/r, 
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tCov/  it. 
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C/5  J 


~To  (_t 
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A 


(4.13) 


The  first  parameter  of  Eq.  (4.18)  may  be  combined  with  the 
second  to  give  the  alternate  parameter  •()£_&_  ) 

Hence  Eq.  (4.18)  may  also  be  written  as  '  L  *"  ' 


fCoy  ~ti 

_T__ 


60s  ~To  L 


A, 


(4.19) 
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ft  ~  ~=SSr  35 

t  A  j  Cps  C'Psj  j  6rs 


Therefore,  it  is  seen  that  with  the  assumption  of  zero  tem- 

“Sr;ienU  ^  th<!  X  ^  y  direCti°"S  *  a  ^e  Juy 

ln  skln>  the  P-irometers  of  Eq.  (4  17)  Ar„  . 

duced,  1^.7)  are  considerably  re- 


If  bota  skins  and  webs  are  thick  the  similarity 
parameters  that  will  eenPrn11v  „  similarity 

and  generally  apply  are  Eqs.  (2.18)  and  (2.21) 

it  is  necessary  to  have  a  model  which  is  a  geometrical'^ 

scaled  replica  of  the  Drof-r>t-vTb«  v  S  ically 

p  L  cne  Prototype.  However,  if  teroperatm-P 

gradients  across  the  thickness  of  f-h«  k 

nicKness  of  the  web  are  negligible,  the 


web  thickness  need  not  be  scaled  in  the  same  manner  as  L>  if 
the  skins  and  webs  are  not  of  the  syne  materials.  This  is 
possible  because  of  the  parameter  of  Eq*  <4*17^ 

which  results  from  the  skin-web  junction  boundary  condition 
and  allows  one  to  use  different  materials  for  the  model  and 

prototype. 

Still  simpler  similarity  relations  can  be  obtained 
for  these  built-up  wing  structures  if  one  is  willing  to  make 
certain  assumptions  regarding  the  heacing-up  process.  During 
the  early  period  of  this  heating  process,  while  the  skin  is 
heating  up,  very  little  heat  is  conducted  into  the  web  and 
the  web  remains  essentially  at  its  initial  temperature.  (This 
assumption  is  further  justified  by  any  joint  resistance  at  the 
skin-web  junction) .  Hence,  the  &eb  will  introduce  no  simi¬ 
larity  parameters  during  this  early  period.  This  situation 
is  applicable  if  the  skin  is  either  thick  or  thin,  and  is  im¬ 
portant  because  during  these  early  times,  the  temperature  dif¬ 
ferences  between  skin  and  web  may  become  large  and  the  result¬ 
ing  thermal  stresses  may  be  appreciable.  The  thermal  simi¬ 
larity  requirements  will  then  be  Eqs.  (4.12),  (4.13)  or  (4,14) 
of  the  immediately  preceding  sections  a)  or  b) . 

At  some  later  time,  it  may  be  assumed  that  the  skins 
have  heated  up  to  close  to  their  steady  state  equilibrium 
temperature  while  the  webs  have  remained  essentially  at  their 
initial  temperature.  Then,  after  this  time,  heat  will  be  con¬ 
ducted  to  the  webs  and  the  skin  temperature  may  be  regarded 
as  a  boundary  condition  at  the  skin-web  junction.  This  as¬ 
sumption  is  again  more  valid  if  there  is  any  thermal  resist¬ 
ance  at  the  skin-web  junction.  The  similarity  parameters  for 
this  situation  at  these  later  times  are  those  required  for 
simulation  of  the  steady  state  temperature  of  the  thin  skin 
and  Eqs,  (4.15)  and  (4.16a)  for  the  web 
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4.3  Stresses  and  Deflections  of  Equivalent  Plates 

The  structural  behavior  of  plates  constitutes  an 
increasingly  important  region  of  structural  analysis  of  air¬ 
craft.  Low  aspect  ratio  wings  are  characteristically  thin 
and  deform  much  like  a  plate.  A  convenient  means  of  repre¬ 
senting  a  wing  of  this  type  is  to  transform  the  actual  struc 
ture  into  an  equivalent  plate  of  varying  bending  and  exten- 
sional  stiffnesses.  The  similarity  requirements  under  the 
combined  effects  of  heating  and  aerodynamic  loading  may  then 
be  derived  from  the  plate  equations. 

The  governing  differential  equations  of  a  heated 
isotropic  plate  are. 


(  9* 
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(4,22) 
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where  p(x,y, t)  is  the  distributed  loading  per  unit  area. 

/  / 

These  two  equations  are  the  Von  Karraan  equations 
for  large  deflections  of  a  plate  with  variable  stiffness 
properties  under  applied  loads  and  thermal  stresses.  The 
quantity  w(x,y ,  t)  is  the  lateral  deflection  of  the  neutral 
surface  from  its  initial  position,  F  is  the  stress  function, 

C  and  D  are  the  extensional  and  bending  stiffnesses  per  unit 
length  respectively,  and  Nj  and  Mp  are  the  midplane  thermal 
force  and  moment  per  unit  length,  respectively.  These  latter 
four  quantities  are  defined  generally  by 


c  = 


/V,  - 
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P  -  \  f  i  Ji 
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Arf-f-  — 


i-v*- 


A 

r 


(A. 23) 


£*  A  Ti  j? 
/-  v 


where  E  is  the  modulus  of  elasticity,  cK  the  coefficient  cf 
thermal  expansion,  V  Foisson  *  s  ratio,  and  A  T  the  temperature 
rise  from  some  initial  value  and  h  is  the  plate  depth.  The 
lateral  loading  p  may  be  due  to  the  combined  effects  of  aero¬ 
dynamic,  inertial,  and  gravity  forces.  As  discussed  in  Section 
A » la), at  supersonic  Mach  numbers  greater  than  about  2.5,  the 
aerodynamic  pressure  loading  p^  may  be  determined  by  piston 
theory,  Eq.  (A, 2),  which  has  the  equivalent  form 


aw 

1% 


>  W\ 
+  tT?r 


(A.2A) 
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The  displacement  w(x,y, t)  will  in  general  include  rigid  body 
displacements  as  well  as  elastic  deformations.  If  the  rigid 
body  displacements  of  vertical  translation  and  pitching  are 
considered,  the  following  equations  of  equilibrium  must  also 
be  satisfied  in  addition  to  Eqs.  (4.21)  and  (4.22) 

-//-  (&  <-i)  J*Jt  *  fj  f*  = ° 

-jf^™  (jp  +  JJ  f*  =  O 

where  m  is  the  airplane  mass  density  per  unit  area.  Intro¬ 
ducing  the  following  non-dimensional  parameters 

w-  w/u9  t  x  -  x/l  , 

t  =  jt/im  ,  F  *  f/c*Ll  ,  C  =  C/C.  ,  P  *  0/p9 

Nt~Nt/Nt0  j  A?7  =  Mr / W~)e  ,  **  -  sr* /srrio 


(4.25) 


(4.26) 


Eqs.  (4.21)  and  (4.22)  in  non-dimensional  form  become 

"  H*‘A' 

~  f* rt,  +  (C*L)\ dF  +  i’p  +  Z  ?F  ifa  1 
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(4.27) 
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The  resulting  non-dimensional  parameters  are 

yf.ML  ,  fr  +  QrtJ)'  ;  ,  /n^V 

Z.  L  p,  t0 1  w. 

J.  ^  .  ar* 

<T»4.  y  /^r.  >  L  rttp  t  Cj  p ,  Mr,  Nr,  /r* ,  &  ,  V 

Po  'Co  P*  &• 

The  first  parameter  is  the  ratio  of  aerodynamic  forces  to 
bending  stiffness  and  is  less  restrictive  than  the  fw/r.e 
parameter  of  Eq.  (2.35)  since  DQ  can  be  varied  by  changing 
the  structural  geometry.  Also  M  and  y  can  be  different. 

The  second  condition  reflects  the  aerodynamic  thickness  effect 
as  expressed  by  piston  theory  and  replaces  the  requirement 
that  M  must  be  simulated  by  itself.  The  third  parameter  is 
the  unsteady  time  parameter  of  the  external  flow,  the  fourth 


Uo 


(4.29) 
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is  the  ratio  of  inertial  loading  to  bending  stiffness  and 
replaces  the  fj,  parameter  of  Eq.  (2.34),  the  fifth  is 

the  ratio  of  gravity  forces  to  bending  stiffness.  The  para¬ 
meter  U*/L  requires  the  neutral  surface  to  deflect  both  simi¬ 
larly  and  in  scale.  However,  if  the  elastic  deformations  are 
small  such  that  the  non-linear  terms  in  w  of  Eq.  (4.28)  are 
negligible,  the  u,/t  requirement  vanishes.  The  seventh  con¬ 
dition  arises  from  the  buckling  characteristics  of  the  wing, 
the  eight  and  ninth  relate  to  the  thermal  force  and  moment 
and  together  replace  the<^#^*^/^  parameter  of  Eq.  (2.34). 
Finally,  the  double-barred  quantities  are  the  distributions 
of  stiffnesses,  thermal  forces  and  moments,  mass,  and  external 
shape  of  the  wing.  No  additional  parameters  will  be  intro¬ 
duced  by  the  non-dimensional  forms  of  Eq.  (4.25). 

By  suitable  combination  with  /o„  ,  U*'/l  , 

and  the  perfect  gas  relationship  «*=  ff/f>  ,  the  parameter 
/y*,tL  to  may  be  rearranged  into  the  alternate  more  familiar 
mass^density  ratio  !/m  ( Similarly,  the  parameter 
yrrt0^-  may  be  changed  by  combination  with  and 

the  above  ,/*(fmL/rn)  to  yield  (u*/l)(  V7fL) .  Using  these  alter¬ 
nate  forms,  the  similarity  requirements  for  the  equivalent 
plate  aeroelastic  problem,  Eq,  (4.29),  can  be  expressed  as 


iff* ml  ,  (f-hi)MJjo  (  v-t*  ,  /  (P»J?)  (yx ) 
p. —  l  '  l  m  i-^ry  '  ‘.ift/ 

c£»_  ,  C.  L"  t  hTu  /  LlMr.  t  C,  P,  Mr,  %rj  -3,  V 
L  Do  Co  (J,P, 


(4.30) 


The  first  parameter,  in  combination  with  A ?=  ? ***/{*,  may  take 
the  alternate  form  9  L  /mo,  where  q  is  the  free  stream  dynamic 
pressure  .  These  specialized  parameters  are  to  be 

compared  with  the  more  general  ones  of  Eq.  (2.35). 
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It  should  be  noted  that  the  parameter  /fb  is  the 
familiar  Froude  number  and  requires  the  gravitational  force 
to  be  scaled  correctly.  This  Froude  number  serves  mainly  to 
define  the  static  equilibrium  position  of  the  model  under 
its  own  weight  and  to  determine  the  equilibrium  angle  of 
attack  in  flying  models.  Since  these  are  generally  not  too 
significant  in  most  aeroelastic  testing,  the  Froude  number 
requirement  can  be  disregarded  (see  Ref.  22). 

The  parameters  of  Eq.  (4.30)  are  also  the  require¬ 
ments  for  flutter  similarity.  The  third  condition,  > 

serves  only  to  define  the  reference  time  tQ  of  the  problem. 

If  sinusoidal  motion  =  we*  ^  is  assumed  as  in  flutter, 
the  reference  time  tQ  can  be  associated  with  the  period 
and  this  similarity  parameter  becomes  the  familiar  reduced 
frequency  .  The  similarity  requirements  for  flutter 

without  any  in-plane  stresses  (such  as  due  to  axial  loads, 

large  deflections  or  aerodynamic  heating)  reduce  mainly  to 

8  —  «*■ 

the  first  five  parameters  of  Eq.  (4.30)  and  also  ,  A,  v  . 

The  conditions  for  static  aerothermoelastic  simi¬ 
larity  are  determined  by  neglecting  the  parameters  xj**/l  and 
'***- Atll  in  Eq.  (4.29).  Comparison  of  these  resulting  static 
similarity  parameters  with  those  of  the  dynamic  case,  Eq. 
(4.30),  indicates  the  absence  In  Eq.  (4.30)  of  the  u'9/l  para¬ 
meter,  the  density  ratio  parameter  /'**•)  and  the  re¬ 

placing  of  the  Froude  irsaber  parameter  by 

The  similarity  parameters  for  a  vibration  model 
are  obtained  by  putt  ag  the  aerodynamic  loading  f9*  equal  to 
zero,  i.e.,  the  first  three  parameters  of  Eq.  (4.29)  will  be 
absent.  For  this  situation  the  important  parameters  relate 
to  the  inertial  forces  and  elastic  restoring  forces  of  the 
wing  and  simulation  of  the  external  flow,  except  as  it  affects 
the  heating  rate  from  the  boundary  layer,  is  not  necessary. 
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The  ***  /o,U  parameter  now  serves  to  define  the  reference  time 
Lo  Problem.  For  sinusoidal  motion  =  of 

vibration,  this  time  tQ  can  be  associated  with  the  period  '/oD 
and  the  similarity  parameters  can  be  rewritten  as 

2-  i  z  , 

'rr’°  4  <*>  ,,  Uc  a?  t  CoL  /Vr0  L  Mr0 

Oo  f  L  *  '  T7  '  -T77p7 

*  =  -  ^  (4.31) 

C  )  D  ,  Mr,  A/  f  ,  syv> 

i 

where  the  parameter  /o  /f  is  the  ratio  of  inertial  forces 
to  gravitational  forces.  The  model,  under  the  conditions  of 
Eq.  (4.31),  will  have  the  sane  non-dimensional  mode  shapes 
at  the  same  non-dimensional  frequencies  L^/p*  as  the  pro¬ 
totype  wing. 


The  stability  of  a  thin  skin  panel  subjected  to 
aerodynamic  heating  and  aerodynamic  pressures  on  one  side  is 
a  problem  of  current  interest.  The  similarity  parameters 
for  the  study  of  the  self  excited  oscillations  and  buckling 
behavior  of  the  panel  may  be  derived  from  section  4.3  by 
properly  modifying  Eqs.  (4.27)  and  (4.23). 


Consider  a  thin  isotropic  plate  of  constant 
ness  b  with  a  uniform  temperature  distribution  across 
thickness.  The  quantities  C,  D,  U  m  Qf  Eq.  (4.23) 
for  this  plate 


thick- 

the 

become 


C  -  *r  $  j  o  _  E  $  nt  =  E <*  AT /  ^  _ 

'  ~"tT  ~  '  7 


(4 • 32) 


The  loading  term  in  Eq.  (4.27)  will  be  modified  by  noting 
that  the  aerodynamic  thickness  term /vMyi, and  the  gravitational 
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loading  are  zero*  The  non-dimensional  parameters  for  the 
large  deflection  behavior  of  the  panel  may  then  be  taken  from 
Eq.  (4.30)  and  reduce  to 


^r(j-)3  'Hr 


Lf/kk),  it,  T),,  V 

17  T 


(4 .  33 ) 


It  should  be  noted  again  that  the  first  parameter  may  have 
the  alternate  form ( $/me6 )  (L/s)  and  the  parameter  serves 

only  to  define  the  reference  time  t  of  the  flow.  The  bound¬ 
ary  conditions  along  the  edges  of  the  plate  will  introduce 
no  additional  parameters  to  Eq.  (4.33).  The  quantities  C, 

D,  N,j,  are  automatically  the  same  for  model  and  prototype  if 
the  same  materials  and  same  temperatures  are  used  for  both. 

It  has  been  shown  in  Ref.  23  by  Hedgepeth  that  for 
cases  where  M  is  greater  than  1.6  and  the  mass  density  ratio 

greater  than  10  (satisfied  by  most  conventional 
panels),  the  aerodynamic  damping  term  ^•/pmL  may  be  neglected. 
This  permits  a  reduction  in  the  parameters  of  Eq.  (4.33)  be¬ 
cause  for  this  situation  the  parameter  (l/tn)ym^A,r»)  drops  out, 
and  the  parameter  ,  is  replaced  by  S'-i,  which  will 

serve  to  determine  the  time  scale. 

For  the  study  of  the  flutter  behavior  of  a  panel 
which  has  already  buckled,  e.g.,  because  of  thermal  stresses, 
it  is  necessary  to  simulate  the  large  deflections  of  the 
plate.  However,  small  deflection  theory  may  be  adequate  in 
predicting  the  onset  of  flutter  of  an  unbuckled  plate.  The 
assumption  of  small  deflections  will  uncouple  Eqs.  (4.27) 
and  (4.28)  and  the  stress  function  F  may  be  solved  for  from 
Eq.  (4.28)  in  the  form 

F -  9)  '■  (4.34) 

Cc 


where  g(x, y)  is  a  function  of  x  and  y  and  depends  on  the 
boundary  conditions.  Substitution  of  F  into  Eq,  (4.27)  re¬ 
sults  in  the  parameter  (a  d  7^  being  replaced  by  (L/s)  . 

Therefore,  the  small  deflection  assumption  removes  the  re¬ 
quirement  that  ^r)c  and  l/$  must  £,e  simulated  independ¬ 
ently  and  replaces  them  with  the  single  requirement  (L/f)(*AT\0 . 

The  similarity  requirements  for  the  buckling  of  a 
panel  by  thermal  stresses  in  the  absence  of  aerodynamic  and 
inertial  loads  may  be  taken  from  Eq.  (4.3;’)  by  neglecting 
the  first  three  parameters.  For  this  static  problem  and  for 
large  deflections  the  similarity  requirements  are 

Uo  L  U  A  t) 

L  $  ( s .  3  5 ) 

If  in  addition  there  is  a  side  loading  ,  such 
as  a  constant  internal  pressure,  the  following  parameter  must 
also  be  simulated  in  addition  to  Eq.  (4.35) 

■f* L  _  ft 

~~cT~  'r/  < 4-->6 ) 

For  this  buckling  case,  small  deflection  theory  may  be  ade¬ 
quate  in  predicting  the  onset  of  buckling,  and  the  parameters 
L/f  and  (a  A  t)0  need  not  be  simulated  independently  but  in  the 
combined  parameter  (t* / , 

jiiL, -Application  to  Low  Aspect -Ratio  Multi-Web  Wing 

The  parameters  derived  for  the  specialized  aero- 
thermoelastic  situations  can  also  be  applied  to  the  low  as¬ 
pect-ratio  wing  of  Fig.  6.  For  simplicity,  assume  the  ribs 
and  webs  are  equally  spaced  a  distance  S  apart,  the  wing  is 
of  uniform  material  throughout,  and  the  skin  and  web  thick- 
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nesses  are  constant.  An  investigation  of  similarity  under 
such  conditions  has  been  considered  by  L.  Ting  in  Ref.  5, 
and  serves  as  an  illustration  of  what  can  be  done  with  spe¬ 
cialized  similarity  cases. 

For  this  low  aspect  ratio  wing  under  the  above 
assumptions,  the  quantities  defined  in  Eq.  (4.23)  become 
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Assume  now  that  the  same  materials  and  same  gas 
are  used  for  both  model  and  prototype,  and  that  the  same  tem¬ 
peratures  exist  throughout  the  model  and  prototype  at  corre- 

if 

spending  instants  of  time.  Letting  CQ,  D0,  denote 

the  quantities  of  Eq.  (4.37)  evaluated  at  some  reference 
cross-sectional  thickness  A 0  >  the  similarity  conditions 

for  dynamic  aerothermoelastic  phenomena,  Eq.  (4.30),  reduce  to 
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These  assumptions  on  the  temperature  distribution  will  be 
re-examined  later. 
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Utilizing  the  1st,  7th ,  and  3th  conditions,  the 
following  relations  between  the  geometry  of  the  model  and 
prototype  can  be  arrived  at, 

(Mq)m  —  (jle)f, 

(fs)«  -  ^y-  (fS)r  I 

(f»)p 

($J)m  -  n  (&*  (Sw)p 
(fs)r 

where  n  is  the  scale  ratio  l'/l 

n  . 

It  can  be  seen  from  the  above  relations,  that  if 
one  is  willing  to  distort  the  model  construction  in  the 
scaling  down  process,  a  pressure  ratio  /(f^p  other 

than  unity  can  be  accomodate i. 

Turning  now  to  the  general  aerodynamic  similarity. 
Section  2,  it  was  seen  desirable  to  reproduce  the  Reynolds 
number  Re  as  well  as  the  Mach  number  M  in  the  flow.  With 
Reynolds  number  and  Mach  number  simulated,  the  model  free 
stream  pressure  and  density  are  scaled  as  follows 


/  (Poo)/*) 


(4 . 40) 


This  pressure  ratio,  together  with  an  additional 
similarity  parameter  J S  introduced  to  provide  for  similar 
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local  buckling  characteristics  of  the  upper  and  lower  skin 
surfaces  (see  Section  4.4,  and  also  Ref.  5),  reduce  the 
relations  between  model  and  prototype  geometry,  Eqs.  (4.39) 


(Jh>)m  —  (&p)p 

(Js)#i  ■=- 

(J*/)*!  =■  >7  P 
Sm  = 


(4.41) 


A  model  so  constructed  would  satisfy  the  aerody¬ 
namic  similarity  requirements  of  Re  and  M  and  also  the  con¬ 
ditions  of  Eq.  (4.30).  A  look  at  the  remaining  similarity 
conditions  is  in  order. 

The  mass  per  unit  area  of  the  structure  shown  in 

Fig.  6  is 


/W  = 


4  S*1 

Sh  J 


(4.42) 


where  PR  is  the  density  of  the  material.  Substituting  this 
into  the  mass  density  ratio  f* 4/W,  condition  gives 


-r 


A.f*  *| 
5  f$  J 


(4.43) 


which  will  be  the  same  for  model  and  prototype  if  (P§)m  =*  (P»)Rt 
i.e.,  similar  materials  for  both.  This  insures  similar  dy¬ 
namic  aeroelasticity  properties. 

The  ^*°Jl  parameter,  as  mentioned  previously,  serves 
to  define  the  dimensionless  time  tQ  associated  with  dynamic 
problems  and  does  not  present  a  problem  in  simulation  here. 


** 


The  (U  parameter  serves  to  define  the  reference 
deflection  uc  and  indicates  that  deflections  will  be  in  the 
same  scale  as  the  model  scale  ratio  rt. 


The  h/ho  parameter  requires  that  the  thickness  dis¬ 
tribution  of  model  and  prototjrpe  be  similar. 

The  Froude  number  parameter  V / 'f  L  will  not  be 
satisfied  under  the  foregoing  conditions  unless  the  scale 
ratio  n  is  equal  to  1.  However,  as  mentioned  previously, 
the  requirement  for  Froude  number  is  not  too  significant  in 
most  aeroelastic  testing  and  can  be  usually  disregarded. 

The  question  of  the  heat  conduction  similarity 
must  now  also  be  considered  in  connection  with  the  original 
assumption  of  the  same  temperatures  existing  throughout  both 
model  and  prototype  at  corresponding  instants  of  time.  Be¬ 
cause  the  model's  dimensions  have  been  distorted  from  those 
of  the  prototype  according  to  Eqs.  (4.41),  no  unique  corre¬ 
spondence  in  temperature  distribution  exists  at  all  periods 
of  time.  However,  under  the  simplifying  assumptions  in  the 
character  of  the  heat  i.low  as  discussed  at  the  end  of  Section 
4.2c),  a  rough  correspondence  in  the  temperature  distributions 
can  be  determined  for  either  very  early  times  or  very  i  *te 
times.  Thus,  at  very  early  times,  it  may  sometimes  be  ade¬ 
quate  to  assume  that  while  the  thin  skin  is  heating  up,  there 
is  no  heat  conduction  into  the  web  and  that  the  web  remains 
essentially  at  zero  initial  temperature  (♦"his  assumption  is 
further  justified  by  any  joint  resistance  at  the  skin-web 
intersection) .  Also,  it  can  be  assumed  that  there  is  negli¬ 
gible  heat  conduction  along  the  skin.  It  can  then  be  shown 
that  the  reference  time  for  this  distorted  model  at  early 
times  is  related  to  that  of  the  prototype  by,  (see  Eq.  (4.13)) 
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This  early  time  assumption  may  be  sufficient  to  cover  the 
case  where  the  maximum  thermal  stresses  are  developed. 

Alternatively,  for  later  times,  it  might  be  assumed 
that  the  thin  skins  have  heated  up  to  the  steady  state  equili¬ 
brium  temperature  value,  and  that  the  heat  is  then  conducted 
into  the  web  according  to  Eq.  (4.15)  subject  to  the  boundary 
condition  of  constant  temperature  at  the  web-skin  junction. 

The  corresponding  reference  times  at  these  later  times  are 
then  related  by  (see  Eq,  (4.20)) 


(  4.45) 


Section  4.5  has  served  as  an  illustration  of  what 
can  be  done  in  the  way  of  satisfying  specialized  similarity 
laws  for  an  equivalent  plate  rather  than  the  most  general 
ones  developed  in  Section  2.  By  distorting  the  geometry  of 
the  model  construction,  it  has  been  possible  to  satisfy  both 
the  M  and  Re  criteria.  However,  the  heat:  conduction  simi¬ 
larity  is  only  approximate  and  contingent  on  certain  assump¬ 
tions  in  the  heat  flow  behavior.  In  view  of  these  latter 
approximations,  it  will  be  shown  later  in  Section  5,  that  it 
is  probably  much  mors  expedient  to  relax  the  Re  requirement 
itself,  rather  than  to  go  to  all  the  trouble  of  geometrically 
distorting  the  model  as  described  here. 
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5 ,  REYNOLDS  NUMBER  RELAXATION 


The  previous  sections  have  discussed  the  problem 
of  satisfying  the  similarity  conditions  by  considering  the 
use  of  different  materials  and  gases,  and  by  determining 
less  restrictive  parameters  by  looking  at  specialized  situa¬ 
tions.  As  mentioned  in  Section  2.4,  a  third  means,  that  of 
relaxing  one  of  the  conflicting  requirements  of  M,  Re,  or 
/Sv  is  also  possible.  The  f* /£■&  requirement  is  essential 
to  any  aeroelastic  type  model  test,  although  it  can  be  made 
somewhat  less  restrictive  by  consideration  of  specialized 
geometries  such  as  equivalent  plates;  see  Eqs.  (4.30).  The 
Mach  number  M  requirement  seems  fundamental  to  most  compres¬ 
sible  flow  studies;  however,  it  too  can  be  made  somewhat  less 
restrictive  by  combination  with  other  parameters  for  certain 
type  flows  as  indicated  in  Section  4.1a)  and  by  Eq.  (4.30). 

This  relaxing  of  M  has  been  suggested  in  Ref.  10  for  certain 
tests.  The  relaxation  of  the  difficult  Reynolds  number  Re 
requirement,  is  also  possible  and  appears  the  most  fruitful 
in  attaining  large  scale  ratios.  As  previously  mentioned 
this  is  often  the  practice  in  ''cold”  aeroelastic  testing, 
provided  a  certain  minimum  Reynolds  number  is  attained.  In 
this  "cold**  case,  simulation  of  the  Reynolds  number  Is  pri¬ 
marily  for  the  purpose  of  simulating  air  loads.  In  the  pres¬ 
ent  problem,  the  Reynolds  number  also  influences  the  rate  of 
heat  transfer  from  the  boundary  layer  into  the  structure  and 
thus  serves  generally  to  define  the  characteristic  time  of 
the  heat  flow  process.  Eq.  (4.7)  shows  how  the  rate  of  heat 
transfer  may  be  approximately  related  to  the  Reynolds  number. 
Thus  the  importance  of  simulating  the  Reynolds  number  may  be 
found  by  investigating  the  effect  of  the  heat  transfer  coef¬ 
ficient  A  »  ft/faw-Tw)  on  particular  aerothermoelastic  phenomena. 

An  important  effect  of  aerodynamic  heating  on  aero- 
elastic  type  models  are  the  resulting  thermal  stresses  in  the 
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structure  which  can  subsequently  cause  reduction  in  the  tor¬ 
sional  rigidity  of  the  wing.  This  loss  in  torsional  rigidity 
results  from  thermally  induced  compressive  axial  stresses  in 
the  leading  and  trailing  edges  and  may  affect  the  structural 
integrity  of  the  wing,  particularly  from  flutter  considera¬ 
tions.  The  problem  of  simulating  the  temperatures,  thermal 
stresses  and  resulting  loss  in  torsional  rigidity  due  to 
these  stresses,  without  simulating  the  Reynolds  number,  will 
now  be  examined  for  some  particular  wing  cross  sections. 


Consider  first  a  two-dimensional  solid  double  wedge 
airfoil  accelerated  instantaneously  from  M0  to  a  final  Mach 
number  Mf  at  constant  altitude.  This  instantaneous  accelera¬ 
tion  results  in  the  greatest  reduction  in  torsional  rigidity 
and  also  is  the  type  of  test,  envisioned  for  aerothermoelastic 
testing.  This  problem  has  been  studied  by  Budiansky  and 
Mayers  in  Ref.  24  which  derives  the  following  expression  for 
the  effective  torsional  stiffness  of  the  airfoil 


Cr  Jo 


/ 


(5.1) 


where  GJ0  is  the  initial  St.  Venant  torsional  stiffness,  <Ty 
the  axial  stress,  negative  if  compressive,  r  the  distance 
from  the  center  of  twist,  and  the  integral  is  to  be  taken 
over  the  cross  sectional  area  A.  The  temperature  in  the  wing 
was  determined  by  assuming  one -dimensional  heat  flow  through 
the  wing  thickness,  no  conduction  in  the  chordwise  direction 
and  constant  material  properties.  At  any  point  along  the 
wing  the  temperature  T  is  given  by  Ref.  24  as 
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(0) 

where  Taw  is  the  adiabatic  wall  temperature  corresponding 
to  Mq,  =  -  7*^  ,  <To  is  the  maximum  thickness  and  A  is 

the  time  parameter  defined  by 

y\  =  z  . 

A  <C 


The  quantities  P g  and  Cpg  are  the  material  density  and  speci¬ 
fic  heat,  respectively  and  t  is  the  physical  time.  The  heat 
transfer  coefficient  h(x)  may  vary  along  the  chord  but  is  not 
a  function  of  time  since  altitude  and  Mach  number  are  constant. 
It  should  be  noted  that  the  parameter  A  corresponds  to  the 
parameter  4* $0L  of  Eq.  (4.13)  for  temperature  similarity 
in  a  thermally  thin  plate.  The  in-plane  thermal  stresses, 
derived  on  the  basis  of  elementary  theory,  have  the  form 

ff  =  Eo(  A  Taw  f  (  >  )  ir)  (5 .  '• ) 


where  x  is  measured  along  the  chord  c.  With  the  use  of  Eqs. 
(5.2)  and  (5.4).  the  effective  torsional  stiffness,  Eq.  (5.1), 
has  the  form 


dr  f-f 

GJ» 


'  -  4rT  An) 


(5.5) 


The  correspondence  is  more  apparent  upon  the  substitution 
of  h(x)  from  Eq.  (4.7)  into  Eq.  (5.3). 


where  A  C'x)  is  a  function  of  ^  and  the  airfoil  geometry. 
Thus,  for  a  given  material,  geometry,  and  initial  and  final 
Mach  numbers,  the  temperatures,  stresses  and  effective  tor¬ 
sional  stiffness  of  the  so? id  double-wedge  airfoil  are  a 
function  of  the  time  parameter  ^  only.  Non-dimensional  tem¬ 
peratures,  stresses,  and  Eq.  (5.5)  will  be  the  same  for  model 
and  prototype  at  corresponding  values  of  ^  .  Since  h  appears 
in  combination  witi  the  time  t  in  the  expression  for  A  ,  an 
incorrect  Reynolds  number,  i.e.,  incorrect  h,  will  merely 
affect  the  physical  time  at  which  the  experimental  results 
occur.  Thus  if  the  model  Reynolds  number  is  lower  than  the 
prototype  value,  similarity  in  T,  &  y,  and  GJr(f/6T0  will  take 
place  at  a  later  time  t  for  the  model  than  for  the  prototype. 
The  discussion  so  far  has  been  restricted  to  an  instantaneous 
acceleration  from  M0  to  M^.  Ref.  24  also  shows  that  for  this 
airfoil  the  maximum,  loss  in  torsional  stiffness,  the  most 
serious  effect  of  the  aerodynamic  heating,  is  relatively  un¬ 
affected  by  moderately  rapid  accelerations.  The  effect  of 
Reynolds  number  relaxation  then  will  be  approximately  the 
same  for  moderately  rapid  accelerations  and  for  instantaneous 
accelerations . 

Consider  next  the  three  cell  built-up  wing  section 
of  Fig.  7.  Assume  the  skin  temperature  to  be  given  by  Eq. 
(5.2)  and  the  web  temperature  at  any  instant  to  be  given  by 


7  w(r)  ? 


C>.6) 


r_  £ 

where  >  0&/z)  •  This  assumes  the  temperature  of  the  skin  is 

unaffected  by  conduction  to  the  web  and  is  valid,  as  mentioned 
in  Section  4.3,  for  skin-web  joints  of  high  thermal  resistance. 
Applying  the  same  assumptions  as  for  the  double-wedge  profile, 
with  h  constant  and  for  n  positive,  the  stresses  and  the 
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effective  torsional  stiffness  for  this  cross  section  are 


£  = 

AT**  ( f~c  )  f  ) 

7  K  n  vy  J 

(5.7) 

= 

(5. a) 

GTtff 

-s  /  -4-  4  T A*  (p 
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where 


Again,  as  for  the  double-wedge,  the  Reynolds  number  appears 
through  the  time  parameter  ^  in  the  temperatures,  stresses, 
and  /crT*  expressions.  The  most  serious  situation  for  the 
built-up  section  regarding  thermal  stresses  occurs  when  the 
webs  are  at  their  initial  temperature  while  the  skin  is  being 
rapidly  heated.  This  corresponds  to  the  case  where  n 
then  reduces  to 


2 


( *■ 


(3 .  ]  1  ) 


When  /$  —  I  ,  in  Eq.  (5,»1!)  is  never  negative,  therefore 
this  built-up  section  will  experience  no  reduction  in  / fC  fa  , 
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However*  loss  in  torsional  stiffness  is  quite  likely  if 
the  skin  panels  buckle  locally  under  the  action  of  thermal 
stresses  induced  by  large  temperature  differences  between 
the  skins  and  webs. 

The  loss  in  torsional  stiffness  of  a  multiweb  air¬ 
foil  with  a  double-wedge  profile  was  studied  in  Ref.  25*  in 
which  heat  was  conducted  to  the  webs  from  the  skin.  The  some 
assumptions  regarding  flight  history  and  skin  temperature 
were  made  as  for  the  preceding  two  examples  and  again  the 
stresses  and  G ****/£!•  were  found  to  depend  only  on  ^  and  geom¬ 
etry.  Fig.  8  shows  210)  calculated  for  this  section  for 
various  ratios  of  the  web  area  to  skin  area  and  a  heat  trans¬ 
fer-coefficient  constant  along  the  chord.  The  results  show 
that  the  loss  in  stiffness  is  less  for  the  skin-web  section 
than  for  the  solid  double-wedge.  In  the  limiting  case  where 
there  are  no  webs*  Aw/tis  =  0,  the  section  becomes  a  hollow 
double-wedge  and  there  is  no  loss  in  stiffness.  This  is  a 
consequence  of  the  constant  heat-transfer  coefficient  assump¬ 
tion.  Fig.  9  shows  A  ( x)  for  a  hollow  section  with  a  vari¬ 
able  h  across  the  chord*  for  laminar  and  turbulent  flow. 
Though  the  reduction  in  stiffness  is  greater  when  h  is  vari¬ 
able  rather  than  constant  across  the  chord*  it  is  still  sub¬ 
stantially  less  than  the  reduction  in  stiffness  for  a  solid 
double-wedge. 

In  summary*  for  the  specific  cases  considered  here, 
under  the  assumptions  of  rapid  accelerations  at  constant 
altitude,  one-dimensional  heat  flow,  elementary  structural 
theory*  and  constant  material  properties*  the  importance  of 
the  Reynolds  number  with  respect  to  temperatures*  thermal 
stresses*  and  loss  in  torsional  stiffness  resulting  from 
these  stresses*  is  manifested  in  the  time  parameter  ~X -  %  • 

A  distortion  in  the  Reynolds  number  only  affects  the  time  t 
at  which  the  experimental  results  occur.  Thermal  similarity 
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thus  is  not  affected. 

The  above  results  regarding  the  influence  of 
Reynolds  number  Re  on  specific  models  can  be  generalized  to 
other  more  arbitrary  configurations*  For  general  thin  skin 
built-up  type  wings,  the  assumption  of  the  skin  heating  up 
one -dimensionally  with  not  much  heat  going  into  the  webs 
during  early  times  (see  discussion  in  Sections  4.2c)  and 
4.5),  can  be  used  to  confine  the  incorrectly  scaled  Reynolds 
numbers  roughly  to  a  redefinition  of  the  reference  time  tQ 
according  to  the  parameter. 


J  to 

Pa 


h 


332  /(£)  (Pr)  4,4 

P,  Cp  f„  L 


(laminar) 

(5.12a) 


f/i-  //, 

=  ^  fcW  (Pr)  Ao 

Pa  Cfs  Sos  L* 


(turbulent) 

(5.12b) 


The  aeroelastic  properties  of  geometrically  scaled,  aerody- 
namically  heated  models  would  then  be  similar  at  correspond¬ 
ing  times  defined  by  Eq.  (5.12),  during  the  early  period  of 
the  heating  up  process.  It  is  expected  that  the  maximum 
thermal  stress  effects  will  be  encountered  during  this  period. 

The  validity  of  all  the  previous  similarity  laws 
and  their  relaxations  should  be  checked  by  constructing  sever¬ 
al  models,  each  at  a  different  scale,  and  comparing  their  non- 
dimensional  responses  under  identical  operating  tunnel  con¬ 
ditions. 
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6.  CONCLUSIONS 

The  present  report  has  reviewed  some  of  the  basic 
similarity  parameters  for  the  general  aerothermoelastic  prob¬ 
lem.  The  investigation  was  limited  to  moderate  temperatures 
(below  about  1000°F  -  i.e.,  M<3.5)  and  considered  only  elas¬ 
tic  behavior  of  the  material. 

The  similarity  parameters  for  general  aerothermo¬ 
elastic  phenomena  in  the  above  range  were  developed  in 
Section  2,  and  are  given  by  Eqs.  (2.35).  Certain  inherent 
conflicts  are  present  which  permit  similitude  only  for  a 
scale  ratio  of  1.  The  primary  conflict  is  between  the  M,  Re, 
and  conditions.  For  same  materials  and  same  tempera¬ 

tures,  these  three  conditions  can  be  satisfied  only  by  a 
scale  ratio  of  1.  Other  less  serious  conflicts  involving 
the  characteristic  time  scale  to,  the  radiation  condition 

f  axid  the  Froude  number  u  ffL*  are  also  present, 
but  these  are  generally  relatively  minor. 

To  get  around  these  conflicts,  one  may  use  differ¬ 
ent  materials  at  different  temperatures  in  different  gases. 
This  is  discussed  in  Section  3.  However,  scale  ratios  ^r/Lm 
of  only  about  3  can  be  achieved  in  this  manner  by  considering 
different  materials,  and  scale  ratios  of  about  3  by  consi¬ 
dering  different  gases. 

Another  means  of  getting  around  these  is  to  develop 
less  restrictive  similarity  parameters  by  considering  more 
specialised  situations.  This  is  discussed  in  Section  4. 

When  applied  to  a  plate  type  structure,  less  restrictive  para¬ 
meters  as  given  by  Eq.  (4.30)  can  be  developed.  Also,  the 
consideration  of  heat  conduction  similarity  in  thin  skin  type 
structures  at  early  times  will  lead  to  the  parameters  of  Eq. 
(4*13) •  Conventional  flutter,  vibration,  and  buckling  phe¬ 
nomena  are  also  shown  to  be  characterised  by  these  specialized 
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similarity  parameters*.  In  Section  4,5,  it  is  shown  that 
geometrically  distorted  models  can  be  constructed  fulfilling 
these  less  restrictive  similarity  conditions. 

A  third  approach  to  these  conflicting  similarity- 
requirements  is  to  leave  one  of  the  general  parameters  un¬ 
satisfied.  This  is  discussed  in  Section  5.  Consideration 
of  the  specialized  situations  above  indicate  that  for  geom¬ 
etrically  scaled  aeroelastic  type  models,  it  is  most  feasible 
to  relax  the  Re  requirement.  For  thin  skin  type  models  at 
early  times,  this  incorrectly  scaled  Reynolds  number  resuits 
only  in  a  redefinition  of  the  time  scale  ratio  according  to 
Eq.  (5.12). 

For  certain  non-aeroelastic  type  tests  in  which  the 
aerodynamic  pressure  distribution  is  known  or  can  be  estimated 
in  advance,  the  f*/Eo  requirement  can  be  replaced  by  the 
more  general  Eq.  (2.40),  and  similarity  can  be  achieved  by 
placing  additional  external  surface  loads  on  the  models. 

This  is  discussed  at  the  end  of  Section  2.4. 

The  success  with  which  each  of  the  above  approaches 
may  be  utilized  to  predict  the  performance  of  a  prototype, 
depends  on  the  specific  situation  involved  and  the  unde:  a  end¬ 
ing  of  what  parameters  may  be  of  significance  to  the  specific 
problem.  In  this  connection,  the  use  of  "restricted”  purpose 
models  can  be  of  considerable  help  in  investigating  clearly 
delineated  problems.  The  validity  of  all  the  previous  simi¬ 
larity  laws  and  their  relaxations  should  be  checked  by  con¬ 
structing  several  models,  each  at  a  different  scale  and  com¬ 
paring  their  non-dimensional  responses  under  identical  tunnel 
operating  conditions. 
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Figure  4  Satisfying  the  General  Similarity 
Parameters  M,  K,,  */&,  K./A  ,  «.% 
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Figure  5  Properties  of  Some  Possible  Substitute  G*ee* 
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FIG.  7  THREE  CELL  BUILT-UP  WING  SECTION 
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SOLID  DOUBLE  WEDGE 


(Taken  from  Ref.  2[j ) 


Time  Parameter,  X 
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FIG.  8  TORSIONAL  STIFFNESS  REDUCTION  OF  A 
MULTIWEB  AIRFOIL. 


(Taken  from  Ref.  25) 


CONSTANT 


•6  .0  1.0 

Time  Parameter,  x 


FIG.  9  EFFECT  OF  A  VARIABLE  HEAT  TRANSFER 
COEFFICIENT  ON  TORSIONAL  STIFFNESS 
OF  A  HOLLOW  DOUBLE-WEDGE  AIRFOIL. 
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